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Summary 

Through the empirical correlation of experimen- 
tal data and theoretical analysis, a set of graphs have 
been developed which summarize the inviscid aero- 
dynamics of delta wings at supersonic speeds. The 
zero-lift wave-drag characteristics of delta wings with 
diamond, circular-arc, and NACA modified four- 
digit-series airfoils were determined through the ap- 
plication of a nonlinear computational technique. 
The nonlinear analysis varied substantially from the 
exact linear-theory predictions for all combinations 
of geometry and flow parameters under study. For 
slender wings with highly subsonic leading edges, the 
nonlinear analysis showed that the zero-lift wave- 
drag correlation relationship was maintained; how- 
ever, as the wing geometry became nonslender, the 
flow about the wing became nonlinear and the rela- 
tionships which define the zero-lift wave-drag corre- 
lation parameter were not maintained. 

The aerodynamic characteristics of delta wings at 
lifting conditions have been evaluated for the effects 
of wing leading-edge sweep, leading-edge bluntness, 
and wing thickness and camber and then summarized 
in the form of graphs which may be used to assess the 
aerodynamics in the preliminary design process. Em- 
pirical curves have been developed for the lift-curve 
slope, nonlinear lift effects, maximum lift, longitu- 
dinal stability, and distribution of lift between the 
upper and lower surfaces of a wing. In addition, the 
impact of various airfoil parameters, wing leading- 
edge sweep, and lift coefficient on the drag-due-to-lift 
characteristics has been shown theoretically. 

The various graphs which detail the aerodynamic 
performance of delta wings at both zero-lift and 
lifting conditions were then employed to define a 
preliminary wing design approach in which both the 
low-lift and high-lift design criteria were combined to 
define a feasible design space. 

Introduction 

An extensive survey of the literature reveals that 
supersonic aerodynamic research began in the 1940’s. 
The literature survey showed that during the 1940’s 
and 1950’s most of the experimental studies at super- 
sonic speeds were conducted on simple wings, bod- 
ies, and wing-body geometries as the aerodynami- 
cists attempted to determine the nature of the new 
flow regime. However, with the introduction of linear 
theory in the 1950’s and the emergence of high-speed 
computers in the 1960’s, there was a gradual de- 
cline in the amount of experimental research directed 
toward the study of fundamental flows over simple 
wings and bodies and a large increase in research di- 
rected toward the application of linear theory to the 


design of aerodynamic vehicles. This change in re- 
search emphasis was no doubt highly influenced by 
the challenge of a new frontier — efficient supersonic 
flight. 

In the design of an optimum wing at super- 
sonic speeds, the aerodynamicist may select from 
an infinite array of wing planforms, airfoil profiles, 
and camber and twist distributions to meet the de- 
sired aerodynamic performance requirements. His- 
torically, the aerodynamicist relied heavily on linear 
theory for the selection of all wing parameters. How- 
ever, comparison of linear-theory estimates with ex- 
perimental data for wings has shown an inability of 
the theory to consistently predict the measured re- 
sults (ref. 1). Linearized theory computational tools 
were limited by the constraints inherent to small 
perturbation theory. 

An extensive survey of the literature was con- 
ducted to determine the dominant wing geometric 
characteristics and flow conditions which should be 
used in assessing the supersonic aerodynamics of 
wings. The initial result of this effort was the identifi- 
cation of the delta or triangular wing planform as the 
most likely candidate for future parametric super- 
sonic wing studies due to the extensive experimental 
data base which was available. In addition, the em- 
pirical correlations derived for delta wings could be 
extended to other simple wing planforms, such as 
arrow and diamond wings, through the use of the 
geometric and flow correlation parameters. 

This paper presents the results of a combined ex- 
perimental and theoretical study of the aerodynam- 
ics of delta wings at supersonic speeds along with a 
preliminary delta wing selection and design philoso- 
phy. The supporting information for this study was 
derived from both the application of a nonlinear in- 
viscid computational method (ref. 2) and from pub- 
lished force, pressure, and flow visualization data for 
delta wings. 

Results of this study should provide a better un- 
derstanding of the effect that airfoil profile, wing 
sweep, and Mach number have on the aerodynam- 
ics of delta wings at supersonic speeds and should 
establish the limitations of linear and nonlinear the- 
ories to predict these effects. The results of the study 
are presented in a parametric fashion in an effort to 
consolidate the effects of wing sweep, Mach number, 
and airfoil shape. 

Symbols 

A wing aspect ratio 

b wing span 

Cd drag coefficient, 



A Cd 

change in drag coefficient relative to 
flat wing at zero lift 

m n 

component of Mach number normal 
to wing leading edge, M cos A 

AC d /0C 2 l 

drag-due-to-lift parameter 

m 

(1 + sin 2 a tan 2 A) 1 / 2 
position of airfoil maximum thick- 

Cd,w 

zero-lift wave-drag coefficient 

ness, expressed as a fraction of local 

C d ,w{DX) 

axial integration drag coefficient 


chord 

C D,w( x ) 

sectional drag coefficient 

N Re 

Reynolds number 

Cl 

lift coefficient. ^ 

Q 

R 

dynamic pressure 

leading-edge radius parameter for 

A C L 

nonlinear incremental change in lift 
coefficient with respect to linear lift 
coefficient, Ci \ a= 2 o° - 20 C La 


NACA modified four-digit-series 
airfoil, Leading-edge radius = 
c[1.1019(ri?/6) 2 ] 

Cl, max 

maximum lift coefficient 

r 

position of airfoil maximum 

Cl 0 

lift-curve slope, evaluated at zero 


thickness 


lift, per deg 

S 

wing reference area 

Cm 

pitching-moment coefficient, 
Pitching moment 

S{x) 

local wing cross-sectional area 


qSc 

t 

wing airfoil thickness 

ACm 

nonlinear increment in pitching- 
moment coefficient with respect to 

X 

streamwise coordinate 


linear pitching-moment coefficient, 

%cp 

streamwise center of pressure 


C m \a=20 o - 20^ CL | a= o o 

y 

spanwise coordinate 

Cn 

r<t 

U N 

normal-force coefficient 
wing lower surface normal-force 

Vcp 

spanwise center of pressure on the 
wing semipanel 

coefficient 

VLE 

spanwise position of leading edge 

° N 

wing upper surface normal-force 

z 

vertical coordinate 

coefficient 

a 

angle of attack 

c p 

coefficient of pressure 

a N 

angle of attack normal to wing 

Cp,o 

coefficient of pressure at zero lift 


leading edge, tan -1 

C* u 

minimum wing upper surface 

0 

= yfW - 1 


pressure coefficient 

0C La 

lift-curve-slope parameter, per deg 

Cp,v 

vacuum pressure coefficient, — 2/^M 2 

/3 cot A 

leading-edge sweep parameter 

C 

wing chord 

1 

ratio of specific heats, 1.40 for air 

C 

wing mean geometric chord 

V 

fraction of local wing semispan 

C r 

wing root chord 

Vv 

semispan location of vortex action 

DX 

marching step size used in nonlinear 


line 


analysis 

A 

wing leading-edge sweep angle, deg 

dCm/dCb 

longitudinal stability level evaluated 
at zero lift, per deg 

T 

airfoil thickness parameter, ex- 
pressed as fraction of local chord, 

K 

zero-lift wave-drag correlation 
parameter, Cpy/ /r 2 A 

Subscripts: 

t/c 

LE 

leading edge 

C 

cross flow 

M 

Mach number 

ref 

reference 
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Discussion 

The present study is directed toward understand- 
ing the fundamental aerodynamic performance of 
delta wings through empirical correlations of experi- 
mental data and nonlinear computations. Nonlinear 
analysis is used only when there are insufficient ex- 
perimental data with respect to a particular geo- 
metric or flow parameter under investigation. The 
data base used in this study is from references 3 
through 25 and summarized in table I. 

The geometric and flow parameters under in- 
i vestigation are presented first. The zero-lift drag 
results are then presented followed by the lifting 
characteristics. 

I 

Aerodynamic and Geometric Parameters 

j The identification of an adequate correlation pa- 

rameter is critical to any empirical study. For delta 
wings, linear theory indicates that at supersonic 
| speeds all aerodynamic characteristics are a function 
of the leading-edge flow condition, as described by 
the leading-edge sweep parameter (3 cot A. In addi- 
tion, the aspect ratio of a delta wing is defined as 
being four times the cotangent of the wing leading- 
edge sweep angle. (See fig. 1.) 

Presented in figure 1 is a plot representing the 
range of leading-edge sweep under study, and pre- 
sented in figure 2 is a plot of Mach number against 
I /? cot A for a range of A. Wings with aspect ra- 
tios varying from a minimum of 0.5 to a maximum 
of 4.0 were selected as being representative super- 
! sonic planforms (fig. 1). Also indicated in figure 2 
1 is the sonic leading-edge condition ( (3 cot A = 1.0) 
for delta wings. The data of figure 2 show that in- 
creasing aspect ratio significantly reduces the Mach 
number for a given value of /? cot A. 

[ The zero-lift drag study concentrated on the ef- 
fect of airfoil profile on the supersonic aerodynam- 
ics of delta wings. The study includes a parametric 
| investigation of diamond, circular-arc, and NACA 
' modified four-digit-series airfoils (ref. 26) in which 
the effect of thickness ratio r for all airfoils, maximum 
j thickness position m for the diamond and NACA 
modified four-digit-series airfoils, and leading-edge 
j bluntness R for the NACA modified four-digit-series 
| airfoil is presented. The range of geometric param- 
eters is graphically presented in figures 3, 4, and 5. 

| Thickness ratio varies from 0.02 to 0.10; maximum 
! thickness position, from 0.2 to 0.5; and leading-edge 
\ radius parameter, from 0.0 to 8.0 within the present 
i study. 


Zero-Lift Wave Drag 

Supersonic linear theory suggests that the zero- 
lift wave-drag (Cpy/) of a delta wing varies as the 
airfoil thickness- to-chord ratio squared (r 2 ) and as 
the wing aspect ratio (A) for a given airfoil pro- 
file. In 1957, Bishop and Cane (ref. 27) incor- 
porated the dependence of the zero-lift wave drag 
upon Mach number into the standard linear-theory 
drag curves by combining the Mach number param- 
eter f3 with the wing geometric parameter A and 
then plotting the linear-theory-dependence parame- 
ter {Cd,w /t 2 A) against the leading-edge sweep pa- 
rameter (f3 cot A). This method of presentation rep- 
resents a plot of zero-lift wave drag against Mach 
number for a given wing airfoil shape rather than 
curves showing the effect on wave drag of one or 
more wing geometric parameters as had been done 
previously (ref. 28). 

The primary intent of this section of the paper 
is to establish the dependence of the zero-lift wave 
drag for delta wings on r 2 , A, and M for various air- 
foil profiles. The dependence is established through 
nonlinear aerodynamic analysis with the method of 
reference 2 because of a lack of experimental data. 
Nonlinear solutions are only presented for subsonic 
leading-edge conditions due to a restriction in the 
present full-potential solution technique. All zero- 
lift wave-drag plots are in the format established by 
Bishop and Cane. 

Diamond airfoil . Computations of the zero-lift 
drag of three-dimensional wings with diamond air- 
foils have dominated the theoretical wave-drag stud- 
ies. The simplicity of the geometry lends itself to 
a closed-form solution of the linearized equations. 
A comparison of the nonlinear and linear solutions 
for delta wings with diamond airfoils is presented 
in figure 6. Nonlinear solutions were obtained with 
the method of reference 2, and the linear solutions 
are taken from the closed-form solutions presented 
in reference 27. Attempts to exactly duplicate the 
closed-form linear-theory solutions of reference 27 
with linear-theory computer codes (ref. 29) were un- 
successful because of the inaccuracies inherent to the 
numerical techniques employed to solve the equa- 
tions. The best results received from the linear- 
theory codes had a maximum error, referenced to 
the closed-form solution, of 5 percent. Nonlinear- 
and linear-theory curves are presented for positions 
of maximum airfoil thickness (m) of 0.2, 0.3, 0.4, 
and 0.5. The linear-theory curves represent the zero- 
lift wave drag of all delta wings. The nonlinear curves 
were computed for a 4-percent-thick delta wing of 
aspect ratio 1.0. The delta wing with aspect ratio 
of 1.0 was selected for nonlinear calculations because 
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of the large variation in Mach number which could be 
analyzed for values of (3 cot A from 0.125 to 1.0. A 
review of figure 2 shows that values of /? cot A be- 
tween 0.125 and 1.0 correspond to a Mach number 
range of 1.12 to 4.12 for a wing with A = 1.0 com- 
pared with a Mach number range of 1.01 to 1.67 for 
a wing with A — 3.0. 

Linear-theory zero-lift wave-drag predictions 
(fig. 6) show a large effect of the position of airfoil 
maximum thickness and increasing /? cot A. The re- 
sults indicate that, for values of / 3 cot A less than 0.6, 
a forward shift in the position of airfoil maximum 
thickness from 0.5 to 0.2 results in a reduction in 
drag of 80 percent. For values of /? cot A greater 
than 0.8, a rearward shift in the position of airfoil 
maximum thickness from 0.2 to 0.5 would also pro- 
duce a reduction in drag of 80 percent. For values 
of (3 cot A between 0.6 and 0.8, the relationship be- 
tween the four linear-theory curves varies drastically. 
In this range of /? cot A, the airfoil maximum thick- 
ness lines become supersonic (Mach angle exceeds the 
maximum thickness line sweep angle) for each airfoil 
at different values of /? cot A. When the Mach num- 
ber normal to the airfoil maximum thickness sweep 
line becomes supersonic, a singularity arises in the 
linear solution which produces a large drag rise (first 
peak). The second peak of each linear- theory curve 
corresponds to the drag rise associated with the oc- 
currence of a sonic leading-edge condition; this occurs 
at j3 cot A = 1.0 for delta wings. As Mach number 
is increased beyond the sonic leading-edge condition, 
j3 cot A > 1.0, linear theory predicts a gradual reduc- 
tion in the drag coefficient for all diamond airfoils. 

Nonlinear zero- lift wave-drag estimates for the 
effect of the position of maximum airfoil thickness 
and Mach number are similar to the linear pre- 
dictions; however, the details of the four nonlinear 
curves are significantly different from their linear- 
theory counterparts. Both linear and nonlinear the- 
ories predict a crossover in the four curves; however, 
the linear-predicted crossover occurs over a range of 
j3 cot A from 0.6 to 0.8, and the nonlinear-predicted 
crossover for all curves occurs at (3 cot A = 0.6, indi- 
cating that, at /3 cot A = 0.6, the drag of a delta wing 
with a diamond airfoil is independent of the position 
of airfoil maximum thickness. On either side of this 
crossover point in the drag curves the effect of the po- 
sition of maximum airfoil thickness is similar to that 
predicted by linear theory. In general, the nonlinear 
analysis results show a much smoother variation with 
/? cot A (no peaks or valleys) and lower drag than the 
linear-theory estimates. Nonlinear theory predicts 
lower drag for all conditions except for the curve for 
m = 0.2 for values of (3 cot A between 0.2 and 0.7 


and the curve for m = 0.3 for values of (3 cot A be- 
tween 0.275 and 0.55. 

The large differences between the linear and non- 
linear predictions do not solely establish the existence 
of nonlinear flow over the wings, such as large shocks 
or expansions. A comparison of the two methods 
shows that the linear method tends to predict a sharp 
pressure gradient due to local surface discontinuities, 
whereas the nonlinear method predicts a much more 
gradual change in pressure. For (3 cot A < 0.4, the 
flow calculated with the nonlinear theory is sub- 
critical in the cross-flow plane ( Mq < 1.00), and 
the forces predicted by linear and nonlinear theo- 
ries are in close agreement. At /? cot A > 0.75, the 
two solutions diverge and a comparison of the cal- 
culated flow fields, with the nonlinear code, reveals 
the development of supercritical cross-flow conditions 
(Me > 1.00). However, for values of /? cot A be- 
tween 0.4 and 0.75, the flow about the wing is sub- 
critical yet the differences between the linear and 
nonlinear curves are significant. To provide insight 
into the flow condition in this region, spanwise sur- 
face pressure distributions from the nonlinear anal- 
ysis are shown in figure 7 for a delta wing with 
a 4-percent-thick diamond airfoil and aspect ratio 
of 1.0 at a Mach number of 2.60, (3 cot A = 0.6. Pre- 
dicted pressure results are presented for airfoils with 
m = 0.2 and 0.5 at streamwise positions of 20, 40, 

60, and 80 percent of the root chord. A comparison 
of the pressures for the two airfoils shows a smooth 
and subcritical character to the flow with larger span- 
wise pressure gradients occurring for the airfoil with 
m = 0.2. The forward position of the airfoil max- 
imum thickness creates an effectively blunter airfoil 
(m = 0.2); this results in larger compression pres- 
sures at the leading edge and lower expansion pres- 
sures as the flow expands around the airfoil ridge line 
when compared with the airfoil with m = 0.5. The 
equivalent drag values, at /? cot A = 0.6, which re- 
sult from an integration of the pressure data are due 
to compensating drag characteristics in the stream- 
wise direction. Over the forward portion of the 
wings (0 < x/c r < 0.3), the local drag of the air- 
foil with m = 0.2 is greater than that of the air- 
foil with m = 0.5; between x/c r of 0.3 to 0.7, the 
drag is similar; and over the aft portion of the wing i 
(0.7 < x/c r < 1.0), the drag of the airfoil with j 
m = 0.5 is greater because of the increased slope j 
at the trailing edge. 

Selected nonlinear analysis was performed next 1 
to establish the dependence of the wave drag upon ] 
the thickness-to-chord ratio squared (r 2 ) and aspect ] 
ratio (A) for delta wings with diamond airfoils. Pre- | 
sented in figure 8 is the predicted nonlinear effect for 


4 



a thickness variation from 0.02 to 0.10 and an as- 
pect ratio variation from 0.5 to 3.0. Results are pre- 
sented as a ratio of the nonlinear-predicted zero-lift 
wave drag for each geometry to a reference nonlinear- 
predicted zero-lift wave drag which has been scaled 
by the ratios of the computed to reference wing 
aspect ratios (A/A re f) and thickness-to-chord ratio 

squared ((r/r ref ) 2 ): 

Cd,w{A,t) 

Cd,W (^ref T ref) (^/^ref) C^Aref) 

Thickness effects, presented on the left of figure 8, 
were evaluated for a reference delta wing with A = 
1.0, m = 0.5, and an airfoil with r = 0.04 at /3 cot A 
= 0.25, 0.5, and 0.75. The nonlinear analysis shows 
a reduced drag with increasing thickness compared 
with the value estimated with the zero-lift wave- 
drag correlation parameter ( K ). The data on the 
left of figure 8 show that, at /? cot A = 0.25, the 
linear-theory-dependence relationship is maintained 
within 2 percent; for f3 cot A = 0.5, a 10-percent 
variation is observed; and for /? cot A = 0.75, a 
20-percent variation is found. At f3 cot A = 0.75, the 
predicted cross-flow Mach number is supercritical for 
the 4-percent- thick wing (figs. 7 and 8); an increase 
in thickness would only tend to magnify these effects. 
Note that solutions could not be attained for wings 
with zero thickness (r = 0.00) with the selected 
nonlinear methodology; however, by definition an 
exact agreement between the linear and nonlinear 
solutions would occur. Presented on the right of 
figure 8 is the effect of aspect ratio for a reference 
delta wing with A = 1.0, r = 0.04, and m — 0.2 
at (3 cot A = 0.5 and 0.75. Nonlinear analysis shows 
a nonlinearly increasing drag with increasing aspect 
ratio and a maximum perturbation of 5 percent for 
/? cot A = 0.75 and 10 percent for j3 cot A = 0.5. 

Nonlinear analysis of delta wings with diamond 
airfoils has shown that for f3 cot A less than 0.5 
(subcritical range), the zero-lift wave drag varies with 
the thickness squared and aspect ratio, but for higher 
values of /? cot A (supercritical range), the nonlinear 
results diverge significantly from the linear-theory- 
dependence relationship. Within this supercritical 
range, the nonlinearity of the flow over the wing 
increases significantly and is shown to be dependent 
upon both the wing and airfoil geometries. 

Circular-arc airfoil. The circular-arc profile is 
a unique class of airfoils which is often used in the 
design of supersonic vehicles. The uniqueness of 
the circular-arc airfoil lies in the fact that it has 
only a single design parameter — airfoil thickness. 


Except for the diamond profile, it is probably the 
airfoil most frequently used in fundamental studies 
of wings. Presented in figure 9 is a comparison of 
the axial distributions of area and drag for a series 
of equivalent 4-percent-thick sharp airfoils on a delta 
wing with A = 1.0 at M = 1.80. The three airfoil 
profiles are a diamond with m = 0.5, a circular arc, 
and a sharp modified four-digit series with m = 0.5. 
The diamond and modified four-digit-series airfoils 
are presented in this section of the paper to provide 
a reference point for discussions on the circular-arc 
profile. The axial distributions of area and drag 
show that the circular-arc airfoil and the equivalent 
sharp modified four-digit-series airfoil (m = 0.5, 
R — 0) are similar, and both have a greater volume 
and a different axial drag distribution than those for 
the diamond airfoil. The drag characteristics show 
that the circular-arc airfoil has higher drag than the 
diamond airfoil at the apex and trailing edge; this is 
due to the increased surface slopes in these regions. 

The effect of Mach number on the nonlinear- 
predicted drag characteristics of the circular-arc air- 
foil is presented in figure 10. Also shown in the figure 
are the nonlinear drag characteristics of the equiva- 
lent modified four-digit-series airfoil and the diamond 
airfoil, along with the linear-theory solution for the 
diamond airfoil with m = 0.5. Nonlinear analysis 
shows a smooth variation in drag for (3 cot A be- 
tween 0.125 and 1.0. Comparison of the drag levels 
shows that the diamond airfoil has the lowest drag, 
the modified four-digit-series airfoil has the highest, 
and the circular-arc profile has drag slightly less than 
that for the modified four-digit-series profile. These 
results are consistent with the data presented in 
figure 9. 

Predicted nonlinear effects of airfoil thickness and 
wing aspect ratio for the circular-arc airfoil (fig. 11) 
are similar to those observed for the diamond air- 
foil. The analysis shows a decreasing drag value with 
increasing airfoil thickness and an increasing drag 
value with increasing wing aspect ratio. The data 
also show a greater percentage change in drag for the 
circular-arc airfoil than for the diamond airfoil (fig. 8) 
due to increasing 6 cot A. This increased nonlinear 
variation in the drag is due to an increased amount of 
nonlinear-predicted flow over the wing with increas- 
ing Mach number. The analysis results presented in 
figures 9, 10, and 11 indicate that a different complex 
and nonlinear flow structure exists for the circular- 
arc airfoil than for the diamond airfoil. In addition, 
the characteristics for the circular-arc airfoil were 
shown to be very similar to those for the equivalent 
modified four-digit-series airfoil; this suggests that 
the circular-arc profile can be treated as a subset 
of the modified four-digit-series airfoil family. The 
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analysis clearly shows that the zero-lift drag of sharp 
airfoils with the same values of r and m is not equal; 
the circular arc has 15 to 40 percent higher wave drag 
than the diamond airfoil but only 4 percent less drag 
than the equivalent modified four-digit-series airfoil. 

Modified four-digit-series airfoil. The NACA 
modified four-digit-series airfoil was selected for the 
study of blunt leading-edge airfoils because of the 
flexibility in defining and altering the airfoil pro- 
file. The NACA modified four-digit-series airfoil 
family is defined analytically with the following 
three parameters: airfoil thickness (r), position of 
airfoil maximum thickness (m), and leading-edge 
radius parameter ( R ). The leading-edge radius of the 
modified four-digit series is a function of the leading- 
edge radius parameter ( R ) and the airfoil thickness 
parameter (r) through the relationship: 

Leading-edge radius = c[1.1019(ri2/6) 2 ] 

This airfoil family is used to evaluate the effect of 
leading-edge bluntness in addition to the effects of 
r, m, A, and M on the zero-lift wave drag. Re- 
sults of the nonlinear analysis for the diamond and 
circular-arc airfoils are referenced. To provide insight 
into the geometric character of this class of airfoils, 
the axial distribution of area and drag for various 
values of the parameters m and R is presented in 
figure 12. Nonlinear results are presented for a 4- 
percent-thick delta wing with A = 1.0 at M = 1.41, 
/? cot A = 0.25. The effect of the position of airfoil 
maximum thickness is shown at the left of the fig- 
ure for sharp airfoils. The area distribution shows 
that moving the airfoil maximum thickness forward 
results in a forward shift in wing volume and a de- 
crease in the maximum cross-sectional area. These 
two effects combine to produce a smoother area dis- 
tribution and lower drag. 

A review of the axial distribution of drag for the 
airfoil with m = 0.2 shows longitudinal symmetry, 
whereas the drag of the airfoil with m = 0.5 is 
mostly generated over the final 20 percent of the 
wing length. The integrated drag values (ref. 2), 
which are listed in figure 12, show that the airfoil 
with m = 0.5 has 60 percent higher drag than the 
airfoil with m = 0.2. As Mach number is increased, 
the character of the drag for each airfoil changes 
dramatically as the compressive pressures begin to 
dominate the flow. The apex drag of the airfoil with 
m = 0.2 would be expected to increase significantly 
compared with the airfoil with m = 0.5, and the 
trailing-edge drag of both airfoils would reduce. 

Shown on the right of figure 12 is the effect of 
leading-edge bluntness on the axial distribution of 


area and computed nonlinear drag. Results are pre- 
sented for a 4-percent- thick airfoil with m = 0.5 at 
R = 0, 4.0, and 8.0, which correspond to a leading- 
edge radius, expressed as a fraction of the chord, 
of approximately 0.0, 0.001, and 0.003, respectively. 
The area and drag data clearly show that leading- 
edge bluntness only affects the forward half of the 
wing. Drag data show that increasing airfoil blunt- 
ness produces a localized increase in drag at the wing 
apex followed immediately by a rapid reduction in 
drag. At x/c r = 0.4, the data show a merging of the 
drag data for the three airfoils which then remain 
coincident over the remainder of the geometry. The 
crossover in the drag data of the three airfoils which 
occurs at the apex of the wing produces a canceling 
effect that results in a drag value of approximately 
21 counts (0.0021) for the three airfoils. This inde- 
pendent nature of the drag due to changes in leading- 
edge bluntness is discussed further later. 

The fundamental nonlinear characteristics for the 
modified four-digit-series airfoil are established for 
4-percent- thick sharp airfoils on a wing with A = 1.0. 
The dependence of the drag on the thickness, as- 
pect ratio, and leading-edge bluntness is established 
separately. 

Presented in figure 13 are the nonlinear-theory 
predictions of drag characteristics for sharp modified 
four-digit-series airfoils with maximum thickness lo- 
cations of 0.2, 0.3, 0.4, and 0.5. Results are pre- 
sented for values of /? cot A between 0.125 and 1.0. 
The nonlinear analysis shows trends similar to those 
predicted for the diamond airfoil: smooth variations 
with (3 cot A and a crossover in the drag of all ge- 
ometries over a range of (3 cot A from 0.65 to 0.8. 
The crossover for the modified four-digit-series air- ! 
foil occurs at a value of (3 cot A of approximately 0.7 J 
compared with a value of 0.6 for the diamond air- 1 
foil. The data show that, for values of (3 cot A less J 
than 0.65, moving the maximum airfoil thickness po- 
sition from 0.5 to 0.2 reduces the drag by a maximum ! 
of 50 percent, and for values greater than 0.8, a rear- ' 
ward shift in the maximum airfoil thickness position 
from 0.2 to 0.5 produces a maximum drag reduction 
of 30 percent. 

The zero-lift wave-drag correlation parameter is 
evaluated for the modified four-digit-series airfoil in | 
figure 14. The trends and levels for both the effects i 
of thickness and aspect ratio are similar to the results j 
for the circular-arc airfoil. The nonlinear predictions 
show a maximum variation of 30 percent for thickness 
variations from 0.02 to 0.10 and a 15-percent varia- j 
tion for variations in aspect ratio from 0.5 to 3.0. 

A comparison of these results with those for the 
circular-arc and diamond airfoils shows an increase in 
the nonlinearity of the characteristics with increasing 
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geometric complexity. However, as mentioned pre- 
viously, the discrepancy between the nonlinear pre- 
dictions and the linear-theory-dependence parameter 
estimates is probably due to a combination of an in- 
creased nonlinear flow field and a less accurate reso- 
lution of the geometry within the solution. 

The effect of leading-edge bluntness on the drag 
characteristics of wings is a highly debated and im- 
portant question in the design of wings and complete 
aircraft configurations. In general, previous super- 
sonic wing design studies which were based upon lin- 
ear theory would frequently employ a sharp airfoil 
or small amounts of bluntness on the leading edge of 
the wing in order to minimize the zero-lift wave-drag 
penalty. 

! Presented in figure 15 is the nonlinear-predicted 
effect of airfoil bluntness on the zero-lift wave drag 
of various modified four-digit-series airfoil geome- 
tries. The nonlinear analysis shows that increasing 
1 leading-edge bluntness in combination with increas- 
ing maximum thickness position, increasing thick- 
I ness, or increasing Mach number generally produced 
increasingly larger drag reductions when compared 
with a sharp leading edge, with a maximum drag re- 
duction of 6 percent occurring for /? cot A = 0.75 and 
! R = 4.0. However, the magnitudes of all the varia- 
tions in drag with increasing leading-edge bluntness 
are small compared with the effect of m, M, r, and A 
observed in figures 14 and 15. This analysis clearly 
shows that leading-edge bluntness is the least in- 
fluential geometric parameter on zero-lift wave drag 
I for subsonic leading-edge conditions, despite the in- 
creased nonlinearity in the wing flow field. 

Presented in figure 16 are the spanwise surface 
f pressure coefficient data predicted by nonlinear the- 
ory for two delta wings with A = 1.0 and 4-percent- 
thick modified four-digit-series airfoils with m = 0.5 
and R — 0 and 8.0. The analysis was performed 
at M = 3.16, /? cot A — 0.75, which corresponds 
[ to geometries which have similar drag values (see 
right side of fig. 15) and are relatively indepen- 
dent of maximum thickness position (see fig. 13). 
Predicted pressure distributions are presented for 
x/c r = 0.2, 0.4, 0.6, and 0.8. The data show that, 

[ at x/c r = 0.2, the blunt airfoil experiences higher 
surface pressures over the total span compared with 
► the sharp airfoil. For x/c r values of 0.4, 0.6, and 0.8, 
the blunt airfoil experiences an expansion (low pres- 
; sures) at the leading edge and an overall pressure dis- 
[ tribution which is lower or equal to that of the sharp 
I airfoil. The lower pressures at the leading edge result 
| from the expansion of the flow around the airfoil nose; 

thus, an aerodynamic thrust force is produced which 
I counteracts the drag increase produced by the more 
positive pressures at the apex of the wing. These 


findings also support the data presented in figure 15 
in which the axial distributions of drag for a wing 
with a sharp and a blunt airfoil were shown to cross 
over just aft of the wing apex producing a canceling 
effect of the local drag. 

Comments on zero-lift wave drag . Predicted 
nonlinear zero-lift wave-drag characteristics of delta 
wings with diamond, circular-arc, and NACA modi- 
fied four-digit-series airfoils have been shown to vary 
substantially from the characteristics predicted by 
linear theory. The nonlinear analysis suggests that 
for slender delta wings at small values of /? cot A the 
zero-lift wave-drag correlation relationship is main- 
tained to a higher degree; however, as the wing as- 
pect ratio or airfoil thickness ratio increases, the 
flow about the wing becomes more nonlinear and 
the zero-lift wave-drag correlation parameter is not 
maintained. The delta wing drag curves presented in 
figures 6, 10, and 13 should be adequate for the pre- 
liminary design studies of delta wings (aspect ratio, 
Mach number, and airfoil selection). 

Lifting Characteristics 

The ability to take into account the effect of airfoil 
profile, aspect ratio, and Mach number on the lifting 
characteristics of delta wings is extremely important 
in selecting the proper wing geometry. The predicted 
flow about a delta wing at the zero-lift condition 
has been shown to be a strong function of these 
parameters, and as the wing is taken to a lifting 
condition, the nonlinear effects would be expected to 
increase significantly (ref. 30). The existing design 
philosophy for delta wings, which is based upon 
linear theory, assumes that the lift and drag due to 
lift of a flat wing are only a function of the leading- 
edge sweep angle and Mach number. The lifting 
characteristics are presented in the form of summary 
graphs in which the drag, lift, and pitching-moment 
characteristics are highlighted. 

Lee-side flow characteristics . At zero angle of 
attack, the flow over the upper and lower surfaces 
of the wing is characterized as attached; however, 
as the wing is taken to angle of attack, the lee-side 
flow characteristics can change dramatically (ref. 21). 
In general the lee-side flow over a wing at angle of 
attack can be divided into two categories: attached 
or separated. The development and character of 
each flow type are dependent upon wing planform, 
airfoil profile, and wing camber in addition to the 
free-stream flow conditions. 

The values of M jy and a jy conducive to the ex- 
istence of an attached flow over the lee side of a 
zero-thick delta wing are shown at the top of fig- 
ure 17 (ref. 21). At zero- to low-lift conditions, the 
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flow is subcritical in the cross-flow plane (small cross- 
flow velocity) with the primary flow direction being 
streamwise. As lift is increased, the flow turning an- 
gle about the leading edge becomes greater, thus re- 
sulting in an acceleration of the flow and an increase 
in the inboard flow component. The cross flow will 
eventually turn streamwise as the flow recompresses. 
Further increases in lift increase both the cross-flow 
component of the flow and the resultant recompres- 
sion as shown in figure 17 and will eventually pro- 
duce a cross-flow shock. The occurrence of a cross- 
flow shock system is an indication of the existence of 
nonlinear supercritical-type flow. Further increases 
in lift result in shock-induced separation of the 
boundary layer and the formation of a leading-edge 
bubble. 

The second type of flow which occurs on the lee 
side of a wing is a leading-edge separation system, 
which is characterized by a viscous, rotational mass 
of air that resides inboard of the wing leading edge. 
(See the lower portion of fig. 17.) At low-lift con- 
ditions, a leading-edge bubble develops, and as the 
angle of attack is increased, the bubble lifts off the 
wing surface and evolves into a vortex. Further in- 
creases in angle of attack result in the development 
of a secondary vortex and eventually a vortex with 
shock occurs. For a vortex system, the vortex body 
is connected to the wing surface via the vortex stem 
or feed sheet. The feed sheet is also a viscous flow 
region that emanates from the wing boundary layer 
near the wing leading edge. The shape and position 
of the viscous vortex system are dependent upon the 
flow field external to the vortex which interacts with 
the vortex system until an equilibrium condition is 
established. In the steady-state model, the vortex 
system acts as a physical boundary to the external 
flow field system boundary. The free-stream flow ex- 
pands around the wing leading edge and follows the 
contour of the vortex system undergoing an expan- 
sion followed by a recompression as the flow turns 
about the vortex. 

This flow field external to, yet influenced by 
the viscous vortex system, is termed the “induced 
flow field.” The extent of the induced flow field is 
characterized by a stagnation point or reattachment 
line on the wing upper surface inboard of the vortex 
body. Inboard of this induced flow field is the 
potential flow field, where the flow is attached in a 
streamwise direction. 

The effects of angle of attack, Mach number, and 
leading-edge sweep on the vortex location were de- 
tailed in reference 31. (See fig. 18.) As shown in 
the sketch, the vortex action line is identified as 
the fraction of the local wing semispan r) v . The 
vortex action line is the position at which the 


vortex normal-force vector should be placed to give 
the same wing bending moment as produced by the 
vortex pressures. The data show that the vortex po- 
sition depends mainly on angle of attack and leading- 
edge sweep angle. Typically the vortex moves from 
a location near 80 percent semispan to a location 
near 50 percent semispan as angle of attack changes 
from 4° to 20°. 

The lee-side flow descriptions presented are typ- 
ical for any wing geometry at both subsonic and 
supersonic speeds. However, the magnitudes of 
these individual effects differ substantially with Mach 
number. The most notable of these differences are 
observed in the static pressures measured on the 
wing upper surface. As Mach number increases be- 
yond 1.0, the free-stream static pressure decreases 
rapidly to significantly reduce the aerodynamic im- 
pact of the expansion over the wing lee side. 

Flat delta wings. This section of the paper draws 
upon several previous empirical correlations (refs. 3, 

4, 32, and 33) as well as the extensive data base in 
order to investigate how wing leading-edge sweep, 
Mach number, and angle of attack affect the lift- 
ing characteristics of delta wings. The supersonic 
aerodynamic characteristics of flat uncambered delta 
wings are presented in figures 19 through 37. All 
aerodynamic characteristics are presented in a para- 
metric fashion in order to consolidate a wide range 
of geometric and flow conditions and in an effort to 
provide a set of empirical aerodynamic curves which 
may be used in preliminary supersonic wing design. 

The characteristic lift-curve slope, nonlinear lift 
effects, and maximum lift for flat delta wings are , 
presented in figures 19, 20, 21, and 22. The data 
presented in figure 19 are an extension to the analy- j 
sis of reference 4 in which additional data sets have | 
been evaluated to establish definitive aerodynamic 
characteristics of delta wings. The experimental 
data shown in figure 19 represent a wide range of i 
wing leading-edge sweep, airfoil thickness, and airfoil 
shape; however, the primary distinguishing feature j 
of the data is that of airfoil bluntness, as noted by 
the open and solid symbols. A comparison of the I 
experimental data with linear-theory results shows 
excellent agreement; for values of /? cot A below 0.5 
and between 0.5 and 1.0, the experimental data break ^ 
away from the linear-theory curve and show lift- j 
curve-slope values below those predicted by linear j 
theory. The data of figure 20 show that not until a j 
value of /? cot A of approximately 2.0 do the exper- 
imental data reach the linear-theory level. A close 
examination of the data also shows a division of the 
two experimental data sets with the sharp leading- 
edge data falling below the blunt leading-edge data. 
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The scattering of the two data sets could be due to 
an effective bluntness that may occur for the sharp- 
edge wings which have large thickness or those which 
have the maximum thickness position near the lead- 
ing edge. The change in characteristic lift-curve slope 
with leading-edge bluntness and increasing f3 cot A 
indicates a change in the local wing flow character- 
istics from the linear-theory model. Linear theory 
assumes that all flow disturbances are weak and are 
propagated along Mach lines. For small values of 
/ 3 cot A, the disturbances are weak and the Mach 
cone is an excellent approximation of the bow shock, 
but as the free-stream Mach number is increased (in- 
creasing /? cot A), the differences between the linear- 
theory Mach cone and the finite-strength bow shock 
increase; this indicates the growth of a strong non- 
linear disturbance field. This effect tends to produce 
a change in flow angularity, a loss in energy, and 
a reduction in the Mach number at the wing lead- 
ing edge. Further increases in Mach number, angle 
of attack, or (3 cot A produce greater changes in all 
flow conditions and a reduced lift effectiveness com- 
pared with those for the linear-theory model. The in- 
creased lift-curve slope with increased wing leading- 
edge bluntness is attributed to both the increased 
expansion of the flow about the wing leading edge 
onto the wing upper surface and the increased lower 
surface compression pressures. 

Presented in figure 20 are lift-curve slopes pre- 
dicted by nonlinear theory for delta wings with sharp 
and blunt leading edges. A comparison of the lin- 
ear, nonlinear, and experimental data shows that the 
linear-theory analysis is in very good agreement with 
the experimental data and nonlinear-theory analy- 
sis for values of /3 cot A less than 0.5. For values of 
/? cot A greater than 0.5, the nonlinear analysis and 
experimental data show a lower lift-curve slope than 
that predicted by linear theory. At these large values 
of the leading-edge sweep parameter, the nonlinear- 
theory analysis of wave drag at zero lift (fig. 6) also 
shows an increase in the nonlinear characteristics of 
the flow and a diverging of the two solutions. 

The effect of airfoil bluntness on the lift charac- 
teristics was studied on the modified four-digit-series 
airfoil. Analysis was performed over a range of an- 
gle of attack from 0° to 10° in which the linearity 
of the lift-curve slope was established within 3 per- 
cent. The nonlinear analysis and experimental data 
are in excellent agreement and both compare well 
with linear theory for values of /? cot A less than 0.7. 
(See fig. 20.) For values of /? cot A greater than 0.7, 
leading-edge bluntness results in a lift-curve slope be- 
low that predicted by linear theory but greater than 
that for sharp airfoils. 


The characteristics of lift-curve slope presented ir 
figure 19 were established by evaluating the experi- 
mental data at zero lift. To establish the extent of 
these linear lift characteristics, data were evaluated 
at an angle of attack of 20° and compared with the 
lift coefficients which result from an extrapolation of 
the experimentally determined lift-curve slope at an 
angle of attack of 20° . These results are presented in 
figure 21 as an increment in nonlinear lift as a func- 
tion of the parameter 4/3 cot 2 A. The data clearly 
show that only for combinations of low Mach num- 
ber and high leading-edge sweep (low aspect ratio) 
does nonlinear increasing lift occur (positive values 
of A C L ). 

This observation was initially reported by Brown 
and Michael in 1954 (ref. 34), Squire, Jones, and 
Stanbrook in 1963 (ref. 10), and Squire in 1967 
(ref. 35) and 1980 (ref. 36). However, with this 
present analysis, the bounds of this unique flow 
condition have been established. For values of the 
correlation parameter 4(3 cot 2 A between 0.5 and 1.0, 
nonlinear lift effects are not present; for values of 
4/3 cot 2 A greater than 1.0, nonlinear decreasing lift 
occurs; and for values of 4/3 cot 2 A less than 0.5, 
nonlinear increasing lift is found. 

The maximum lift coefficient for thin delta wings 
plotted against the parameter M/A is presented in 
figure 22. The supersonic aerodynamics of delta 
wings at high angles of attack is presented to ensure 
a complete review of all pertinent aerodynamic char- 
acteristics. Possible applications of the data at these 
characteristics would be to canards or horizontal tails 
of aircraft and fins of missiles. An interesting point to 
note is that, despite the large range in Cl, max which 
was observed in the data (1.1 > Cl, max > 0.7), 
all data showed that the maximum lift occurred at 
an angle of attack between 40° and 50°; this indi- 
cates that Cl, max is dominated by the rotation of 
normal-force vector and not a function of flow con- 
ditions on the wing upper surface as is true at sub- 
sonic speeds. Experimental data at subsonic speeds 
show that wing upper surface flow separation or vor- 
tex breakdown has a large influence on Cl, max- The 
data show that CL, ma x varies from a maximum of 1.1 
at a value of M/A of 0.05 to a minimum value of 0.7 
at a value of M/A of 9.0; this indicates the effect 
of both Mach number and wing sweep. The com- 
bination of high supersonic Mach number and high 
wing sweep which is required to achieve large values 
of the correlation parameter {M/A) corresponds to 
wing flow conditions which would produce reduced 
values of CL, ma x- As an example, experimental data 
show that the contribution to the total lift coefficient 
from the wing windward side is only slightly depen- 
dent of wing sweep and Mach number at maximum 
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lift coefficient. However, the leeward-side contribu- 
tion to the lift coefficient is greatly reduced with in- 
creases in both wing sweep and Mach number; these 
effects combine to produce the characteristics shown 
in figure 22. These individual wing upper and lower 
surface lift characteristics are discussed in more de- 
tail in the discussion of figures 27 through 38. 

The longitudinal stability characteristics of 
uncambered delta wings, presented in figures 23 
through 26, were determined by computing the slope 
across the zero-lift condition. The data have all been 
reduced about the two-thirds root-chord location, the 
linear-theory-predicted moment center. Presented 
in figure 23 is the longitudinal stability parameter 
/3( dCm/dCi ) evaluated at zero lift plotted against 
the leading-edge sweep parameter 3 cot A. Experi- 
mental data are presented for both sharp- and blunt- 
leading-edge wings with various leading-edge sweep 
angles and airfoil shapes over a wide range of Mach 
number. An examination of the data in figure 23 
shows a gradual forward movement in the wing cen- 
ter of pressure when the wing leading-edge flow con- 
dition changes from subsonic to supersonic. For the 
subsonic leading-edge condition, all data show either 
neutral stability or positive stability, whereas for the 
supersonic leading-edge condition, all data show a 
negative stability level. As mentioned previously, 
these data were obtained at the zero-lift condition. 
To evaluate the extent of the linearity of these char- 
acteristics, the zero-lift dCm/dCi, level was extrap- 
olated to an angle of attack of 20° and compared 
with the measured experimental data at that angle 
of attack. Presented in figure 24 are the nonlinear 
increments in pitching moment evaluated at an an- 
gle of attack of 20°. The data show a trend with 
increasing 3 cot A which is opposite that for the sta- 
bility data of figure 23. For fd cot A less than 1.0, 
the data show a positive increment in pitching mo- 
ment, and for fl cot A greater than 1.0, the data show 
negligible increments in pitching moment. The data 
of figure 23 indicate that, for the subsonic leading- 
edge condition at zero lift, the majority of the data 
show that the wing stability level was approximated 
well by linear theory; however, as the wing is taken 
to an angle of attack, the data of figure 24 show 
the development of a nonlinear increment in pitch- 
ing moment for the subsonic leading-edge condition, 
and this suggests the development of nonlinear flow 
over the wing. 

The longitudinal stability characteristics are pre- 
sented in figures 25(a) and 25(b) where the stream- 
wise and spanwise wing centers of pressure on a 
single wing panel are presented for angles of attack 
of 5° and 50°, respectively. The data of figure 25(a) 


show a variation in both the spanwise and stream- 
wise center of pressure with changes in wing aspect 
ratio, whereas the data at an angle of attack of 50° 
(fig. 25(b)) collapse into a single curve. The variation 
in wing center of pressure can be directly related to 
the flow conditions about the wing. At low angles of 
attack, the lee-side flow characteristics for each wing 
vary considerably between attached and separated 
flow; however, as the wings are taken to a = 50°, 
the flow becomes quite similar and is characterized 
by two wing leading-edge vortices. 

The lift and pitching-moment data presented 
have been shown to be dominated by wing leading- 
edge sweep effects. A review of the existing data 
base failed to uncover an adequate parametric set of 
wing data which could be used to isolate the effects 
of geometry and flow conditions on the drag-due-to- 
lift characteristic; therefore, the nonlinear computa- 
tional method of reference 2 was employed for this 
analysis. 

Nonlinear-predicted drag-due-to-lift characteris- 
tics for blunt airfoils are presented in figures 26 
and 27. Also noted in each figure are the linear- 
theory zero-thrust and full-thrust boundaries. The 
nonlinear-theory method selected for this study is 
valid only for conditions of attached flow; as a re- 
sult, nonlinear analysis was limited to blunt airfoils at 
conditions of low-to-moderate lift coefficients (Cl < 
0.3). The effect of airfoil thickness, airfoil maximum 
thickness position, and airfoil bluntness on the drag- 
due-to-lift parameter has been studied for a delta 
wing with aspect ratio 1.0 at a Mach number of 1.41 
(fig. 26). Results of the analysis show that increasing 
thickness, moving the maximum thickness location | 
forward, or increasing airfoil bluntness improves the 
lifting efficiency of the airfoil, with airfoil thickness ' 
providing the largest improvement and bluntness the | 
smallest. Comparison of the results shows that the j 
nonlinear theory predicts a lower drag-due-to-lift pa- 
rameter value than that predicted by linear theory for . 
uncambered wings (no thrust), and for low values of 
the lift coefficient (Cl < 0.1), the nonlinear theory 
predicts drag characteristics below the linear-theory 
optimum (full thrust). The increase in the drag-due- 
to-lift parameter with increasing lift coefficient is a 
result of a loss in the aerodynamic thrust force as the • 
wing rotates through a range of angles of attack. At l 
low-lift conditions, the flow about the wing surface ; 
is characterized by a gradual expansion about the 
nose of the airfoil followed by a smooth recompres- 
sion on the upper surface; this produces significant 
amounts of aerodynamic thrust. As angle of attack is 
increased, the local expansion about the airfoil nose , 
increases significantly and extends over a greater por- . 
tion of the airfoil leeward surface; this results in a 
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smaller percentage of the local expansion pressures 
acting on the nose of the airfoil. This expansion re- 
sults in a reduction in the aerodynamic thrust force of 
the wing. The favorable flow conditions which exist 
about the nose of a blunt airfoil (low pressures) must 
coincide with favorable flow conditions over the aft 
section of the airfoil (high pressures) if an improve- 
ment in the drag-due-to-lift characteristics is to be 
realized. 

The improvement in drag-due-to-lift characteris- 
tics with increasing airfoil bluntness is contradictory 
to the results of previous studies (ref. 37) in which 
complex wing geometries were employed in the study 
of the leading-edge thrust phenomenon. However, as 
discussed, the variation in drag due to lift is not solely 
dependent upon the loading at the leading edge of the 
wing, but it is dependent upon the total wing load- 
ing. The nonlinear analysis presented in figure 26 
was performed at conditions which were previously 
determined to have a subcritical-type flow structure 
over the wing (0 cot A = 0.25); however, the drag- 
due-to-lift analysis showed a highly nonlinear char- 
acteristic. To investigate these effects further, non- 
linear analysis was performed on both a delta wing 
with A = 1.0 and one with A — 2.0 with blunt mod- 
ified four-digit-series airfoils (r = 0.40, m = 0.5, and 
R = 4.0). Analysis was conducted for values of 
0 cot A from 0.25 to 0.8. Results of this analysis are 
presented in figure 27 in which curves of constant 
lift coefficient for both wings are shown. Also pre- 
sented in the figure are the linear-theory zero-thrust 
and linear-theory full-thrust curves of the drag-due- 
to-lift parameter. The nonlinear analyses show an in- 
crease in the drag-due-to-lift parameter with increas- 
ing lift coefficient, decreasing 0 cot A, and decreasing 
aspect ratio. Results for the wing with A = 1.0 at 
Cl =0.1 show a variation in the aerodynamic thrust 
from 70 percent of the full-thrust linear-theory value 
at 0 cot A = 0.25 to 0 percent at 0 cot A = 0.75. The 
drag-due-to-lift parameter for the wing with A = 2.0 
at Cl — 0.1 varied from 114 percent to 100 per- 
cent at 0 cot A = 0.25 to 0.75. Similar character- 
istics to those observed at Cl = 0.1 were also ob- 
served for the two wings at Cl — 0.2. An increase 
in drag due to lift with decreasing wing aspect ra- 
tio would normally be expected for a constant Mach 
number; however, by presenting the data with re- 
spect to 0 cot A, wing leading-edge sweep or aspect 
ratio effects should be factored out of the solution. 
As a result, the analyses suggest that Mach number 
has a large impact on the aerodynamic performance 
of a given wing. As Mach number is reduced, the 
characteristics around the wing change dramatically; 
the upper surface leading-edge suction pressure coef- 
ficients become more negative and the lower surface 


leading-edge pressure reduces slightly to increase the 
aerodynamic thrust force. In addition, the thrust 
force which is produced normal to the wing leading 
edge varies as the cosine of the leading-edge sweep 
angle. As aspect ratio is increased the wing sweep 
decreases to increase the aerodynamic thrust force. 
In general, the character of all the curves presented 
in the figure is a rapid reduction in the drag-due-to- 
lift parameter (increased performance) with increas- 
ing 0 cot A up to a value of approximately 0.6. For 
values of 0 cot A greater than 0.6, the curves tend 
to flatten out. This characteristic seems to be quite 
sensitive to aspect ratio and lift coefficient; however, 
a value of 0 cot A of 0.6 seems to be a representative 
mean value for the analysis conducted. 

The previous analysis reviewed the aerodynamic 
forces and moments of flat delta wings at lifting con- 
ditions in which the lift and pitching-moment char- 
acteristics were shown to be dominated by the flow 
condition at the wing leading edge (0 cot A) and the 
drag to be equally sensitive to all geometric and flow 
parameters. To provide further insight into these ob- 
served characteristics for flat delta wings, wing upper 
and lower surface static-pressure distributions are re- 
viewed. The experimental pressure data are reviewed 
to investigate the effects of wing leading-edge sweep, 
leading-edge bluntness, wing thickness, wing camber, 
Mach number, and angle of attack on the total wing 
loading and the individual upper surface and lower 
surface wing loadings. This study makes use of the 
known existence of conical flow for conical geometries 
at supersonic speeds to extend the use of the limited 
amount of wing surface pressures in an effort to rep- 
resent the total wing loading. However, before this 
analysis can be generally applied, the lower Mach 
number bounds of the supersonic flow regime need 
to be identified in order to determine where the con- 
ical flow assumption may be used. In order to deter- 
mine the lower Mach number bound of the supersonic 
flow regime, subsonic and transonic data need to be 
reviewed to ensure that the characteristics of these 
two flow fields are differentiated from those at su- 
personic speeds. Presented in figure 28 are spanwise 
pressure distributions at various longitudinal ( x/c r ) 
stations for a delta wing with A = 1.0 at Mach num- 
bers from 0.60 to 3.50. The pressure data clearly 
show that nonconical flow conditions exist at both 
M — 0.60 and 0.90. However, between M = 0.90 
and 1.20, the pressure data change from a nonconical 
nature to a very nearly conical flow on both the up- 
per and lower surfaces. At a Mach number of 1.60, 
the spanwise pressure distributions are conical and 
remain conical for all further increases in Mach num- 
ber. Based upon these data, a minimum supersonic 
Mach number is defined as approximately 1.20. An 
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interesting point to note is the dramatic change in 
the upper surface pressures and the comparative in- 
sensitive nature of the lower surface pressures with 
increasing supersonic Mach number. These effects 
are discussed further in the discussion of figure 29. 

Typical effects of Mach number and angle of 
attack on the spanwise pressure distribution for a 
delta wing with A = 1.0 are presented in figures 29 
and 30, respectively. To assist in the discussion of 
the pressure data only a single pressure distribution 
is presented in figures 29 and 30. In addition, the 
data of figures 29 and 30 are used to represent the 
total wing pressure distribution (conical flow). 

If the assumption is made that conical flow exists, 
then the data presented in figure 29 show that, at a 
constant angle of attack of 15°, increasing the Mach 
number from 0.6 results in an ever increasing reduc- 
tion in upper surface normal-force coefficient. A re- 
versal of these upper surface effects is found on the 
lower surface where the lower surface normal-force 
coefficient increases with an increase in Mach num- 
ber as the transition is made to supersonic speeds 
and then remains constant. As Mach number in- 
creases, the increase in lower surface lift coefficient 
does not completely compensate for the reduction in 
upper surface lift, and the result is a reduction in to- 
tal wing lift coefficient with increasing Mach number 
at a constant angle of attack. The combination of a 
reduction in the upper surface lift coefficient with an 
increase in lower surface lift coefficient highlights the 
impact of compressibility and the vacuum pressure 
coefficient at supersonic speeds. 

Presented in figure 30 is the effect of increasing 
angle of attack on the spanwise pressure distribu- 
tion at M = 0.60, 0.90, 1.20, and 1.60 for the delta 
wing with A = 1.0. The pressure data show that at 
M = 0.60 and 0.90, both the upper and lower surface 
lift increase proportionally. The data for M — 1.20 
and for M = 1.60 indicate a shift in the percent- 
age of lift force from the upper surface to the lower 
surface. This effect is quantified in figure 31 which 
depicts the percentage of lift on both the upper and 
lower surfaces for angles of attack of 10° and 20°. 
The graph of figure 31 shows that at subsonic speeds 
(M = 0.60 and 0.90), approximately 70 percent of 
the total lift comes from the upper surface, indepen- 
dent of angle of attack. At M = 1.20 and a = 10°, 
the upper surface carries 70 percent of the lift, but 
at a = 20°, this is reduced to 60 percent. This trend 
is also observed at M = 1.60 where at a = 10° the 
upper surface dominates with 56 percent of the lift 
and at a = 20° it has been reduced to 48 percent 
of the lift. The shifting of lift from the upper to 
the lower surface with increasing angle of attack at 


supersonic speeds is a combination of nonlinearly de- 
creasing upper surface normal force and nonlinearly 
increasing lower surface normal force. These results 
suggest that only at subsonic speeds ( M < 0.90) can 
the lift increment between the linear potential-theory 
solution and experimental data belong solely to up- 
per surface vortex-induced effects. At both transonic 
and supersonic speeds, the lift increment is probably 
due to a combination of both nonlinear lower surface 
and nonlinear upper surface effects. In particular, at 
supersonic speeds, the upper surface vortex-induced 
lift increment reduces with increasing angle of attack, 
and the lower surface compression lift increment in- 
creases with increasing angle of attack. 

At both subsonic and supersonic speeds, the lift- 
curve slope and lifting efficiency of a wing are directly 
related to the aspect ratio of the wing. To determine 
the effect of aspect ratio on the distribution of lift 
between the upper and lower surface of a wing, the 
analysis presented in figures 30 and 31 was repeated 
for a wing with A = 0.5 (A = 82.87°) and one 
with A = 2.0 (A = 63.43°). Presented in figure 32 
are spanwise surface pressure distributions for wings 
with A = 0.5, 1.0, and 2.0 at angles of attack of 
10° and 20° for M — 0.60 and 1.60. Comparison 
of the upper and lower surface pressure plots in 
figure 32 shows that the supersonic data vary in a 
more orderly fashion with increases in aspect ratio. 
The data at both M = 0.60 and 1.60 show an increase 
in lift with increasing aspect ratio, with the upper 
surface dominating at subsonic speeds, but a more 
equal distribution is evident at supersonic speeds. 
These effects are quantified in figure 33 in which the 
percentage of lift on the upper and lower surfaces is 
depicted. The graph shows that Mach number has 
the largest impact on the distribution of lift for all 
wings at angles of attack of 10° and 20°. All wings 
show a shifting of lift from the upper surface to the 
lower surface with increasing angle of attack except 
for the wing with A = 0.5 at M = 1.60, which shows 
a reverse effect. The combination of low aspect ratio 
and low free-stream Mach number creates conditions 
favorable to the development of nonlinear lift effects. 
(See fig. 21.) These conditions combine to give 
a value of the nonlinear lift correlation parameter 
(4(3 cot 2 A) of 0.08, which is a condition at which 
large amounts of nonlinear lift would occur. Based 
upon the data in figure 33 it may be concluded that 
this increment can be attributed to the wing upper 
surface flow field. A comparison of the wings with 
A — 1.0 and 2.0 at M = 1.60 shows a gradual 
increase in lower surface dominance with increasing 
wing aspect ratio at both angles of attack. The 
pressure data of figure 32 show that this transition in 
lift force is primarily due to an increase in the lower 


12 



surface contribution to the normal force and not due 
to upper surface effects. 

A summary of the individual normal-force char- 
acteristics of the upper and lower surfaces for thin 
delta wings is presented in figures 34, 35, and 36. 
Figures 34 and 36 are formulated by integrating 
experimental spanwise pressure distributions to ex- 
tract sectional upper and lower surface normal-force 
coefficients. The sectional normal-force coefficients 
are then used to represent the total wing upper or 
lower surface normal-force coefficient based upon the 
known existence of conical flow for delta wings at su- 
personic speeds. The curves of figures 34, 35, and 36 
cover a range of Mach number from 1.50 to 3.50 and 
leading-edge sweep of 52° to 85°. The results could 
not be produced for Mach numbers below 1.50 be- 
cause of insufficient experimental data. As shown in 
figure 34, when the upper surface normal-force coef- 
ficient Cjy is plotted as a function of the parameter 
aw (3 cot A, the data collapse into a family of constant 
Mach number curves. The large effect of Mach num- 
ber on upper surface normal force is clearly shown; 
for example, an increase in Mach number from 1.50 
to 2.00 (fig. 34(b)) reduces the upper surface lift- 
ing potential by 50 percent at a given value of the 
correlation parameter. The large reduction in up- 
per surface lifting capability with increasing Mach 
number is due to the inability to achieve low val- 
I ues of the upper surface suction pressure coefficient. 
The characteristics of the data of figure 34 also indi- 
cate that, for a given Mach number and leading-edge 
i sweep, an increase in angle of attack results in an 
upper surface normal-force coefficient that increases 
nonlinearly with a decreasing slope. 

Presented in figure 35 is the upper surface mini- 
I mum pressure coefficient C™ min plotted as a function 
of the parameter ajv/3 cot A. These data also collapse 
into a family of constant Mach number curves. Also 
I noted in this figure is the percent of vacuum limit 
• which was attained for that particular Mach num- 
> ber. The data show that the percent of the vacuum 
pressure actually attained is reduced with increasing 
Mach number. At a Mach number of 3.50 (fig. 35(a)), 

| only 75 percent of the vacuum limit pressure co- 
| efficient was reached; however, at a Mach number 
t of 1.50, 97 percent of the vacuum limit pressure 
f coefficient was attained (fig. 35(b)). 
j The variation in the lower surface normal-force 
| coefficient is presented in figure 36. Unlike the data 
j of figures 34 and 35, which showed that the upper 
surface characteristics are a function of leading-edge 
sweep, Mach number, and angle of attack, the lower 
surface characteristics were found to be predomi- 
nantly a function of leading-edge sweep and angle 


of attack. The lower surface normal-force coefficient 
C^ N is plotted as a function of the normal angle of 
attack and the data collapse into a family of 
constant leading-edge sweep curves. Each of these 
curves is comprised for a range of Mach numbers, 
with the only limitation being that the normal Mach 
number M n of all these data be less than 1.00. How- 
ever, the data have a maximum variation with Mach 
number of approximately ±0.03(7/v* The data show 
that the lower surface produces a nonlinearly increas- 
ing normal-force increment with increasing angle of 
attack, and it can be seen that the nonlinearity in- 
creases with an increase in leading-edge sweep. 

An evaluation of the data of figures 34 and 36 sup- 
ports the findings previously observed in figures 21 
and 33, which indicate that nonlinear lift is most 
pronounced for extremely highly swept wings at low 
Mach numbers. For these very low values of the pa- 
rameter /? cot A, the highly nonlinear character of 
the lower surface normal-force coefficient adds to the 
linear character of the upper surface normal-force 
coefficient and produces a total lift force which in- 
creases in a nonlinear sense. Similarly, figures 34 
and 36 can be used to show that the reduction in lift- 
curve slope with an increase in leading-edge sweep is 
primarily a lower surface-dominated effect and the 
increase in lift-curve slope with a decrease in Mach 
number is an upper surface-dominated effect. 

To further extend these findings, the effects as- 
sociated with varying leading-edge bluntness and 
Reynolds number ( Nr 6 ) on the minimum upper sur- 
face pressure are presented in figures 37 and 38, re- 
spectively. The data for C™ min are plotted as a func- 
tion of the parameter ajy/3 cot A for flat wings at 
conditions in which the flow has separated at the 
leading edge. A summary of the effect of leading- 
edge bluntness on C™ min is presented in figure 37(a) 
in which curves are presented for both sharp and 
blunt leading-edge wings at Mach numbers of approx- 
imately 1.50, 1.60, 1.90, and 2.40. The blunt leading- 
edge curves were taken from figures 37(b) and 37(c), 
and the sharp leading-edge data curves were taken 
from figure 35. The data show that leading-edge 
bluntness increases the minimum upper surface pres- 
sure coefficient for a given value of cot A but 
does not affect the minimum value of The in- 

p,min 

creased upper surface suction pressures with leading- 
edge bluntness would be expected to increase the lift- 
curve slope and allow for the development of a given 
level of lift at a lower wing incidence angle, and, thus, 
result in lower drag. These results are directly re- 
flected in the data of figure 19 where the lift-curve 
slope is shown to increase over that for a sharp wing 
with the addition of leading-edge bluntness. 
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The data of figure 38 show that at a constant 
Mach number increasing Nft e for a blunt leading- 
edge wing increases the minimum upper surface pres- 
sure coefficient for a given value of cot A; how- 
ever, these effects decrease with an increase in Mach 
number. 

These data show that at supersonic speeds signif- 
icant amounts of nonlinear flow exist on flat, sharp- 
leading-edge delta wings, and that the nonlinear 
characteristics of delta wings are affected by leading- 
edge bluntness and Reynolds number. These non- 
linear characteristics which occur on both the upper 
and lower surfaces of flat wings are summarized in 
the form of graphs to provide a fundamental under- 
standing of the aerodynamics of delta wings to the 
designer. 

Thickness effects . The figures discussed in the 
previous section, which summarized the global and 
local aerodynamics of flat delta wings, relied heavily 
on the known existence of conical flow for their de- 
velopment. The figures were constructed from data 
sets obtained with zero-thick delta wings (conical ge- 
ometries) or very slender geometries which are near 
conical in nature. In order to study similar aerody- 
namic effects for typical wing thickness variations, 
it would be necessary to have an extensive amount 
of upper surface and lower surface pressure distri- 
butions in order to resolve the nonconical flow condi- 
tions. However, if the thick-wing geometry is conical, 
the flow may be assumed to be conical and a single 
spanwise pressure distribution is all that is required 
on both the wing upper surface and lower surface to 
represent the total wing flow field. 

In an attempt to assess the effect of leading-edge 
angle, a comparison is made between the 30° dia- 
mond wing and the 60° biconvex wing. In an attempt 
to assess thickness effects, a comparison is made be- 
tween the 30° biconvex wing and the 30° diamond 
wing. Presented in figures 39 through 48 are the 
effects of thickness on the total wing aerodynamics 
and the local wing loadings for several conical ge- 
ometries. The thick-wing data were obtained from 
reference 23 in which three delta wings with A = 1.0 
differing in thickness and cross-section shape were 
tested between M = 1.30 and 2.80. As shown in the 
inset sketch of figure 39, the three different thick- 
nesses correspond to a diamond cross section with 
a 30° leading-edge half-angle, measured in the cross- 
flow plane, and two biconvex cross sections with 30° 
and 60° leading-edge half-angles, also measured in 
the cross-flow plane. It is recognized that these cross 
sections do not form typical airfoil shapes; however, 
they probably represent reasonable upper bounds on 
airfoil leading-edge angles and airfoil thicknesses. 


The effect of wing thickness on the wing leading- 
edge vortex position is presented in figure 39. The 
zero-thickness flat-wing curve shown in the figure has 
been taken from the flat-wing data of reference 31. 
The data clearly show an outward movement of the 
vortex with increasing leading-edge slope at all con- 
ditions evaluated. The thick-wing data show that 
the wing leading-edge angle acts to delay the onset 
of flow separation to a higher angle of attack. For a 
constant angle of attack, the data show a weaker vor- 
tex which is located more outboard compared with 
the flat-wing data (fig. 42). A comparison of the 
data for the two wings with a 30° leading-edge sweep 
angle shows a slight shift in the vortex position due to 
thickness. The data indicate that wing leading-edge 
angle and not wing thickness is the dominant mecha- 
nism which controls vortex formation, strength, and 
position. Despite these large effects on the lee-side 
flow characteristics between the three wings, there 
were no noticeable changes on the total wing lifting 
characteristics. (See fig. 40.) For values of /3 cot A 
below 0.5, the lift-curve-slope data show an insensi- 
tivity to thickness; for values of P cot A above 0.5, 
the data show an increase over the flat delta data 
but yet show no variation between the three thick 
wings. The lift-curve-slope data were evaluated be- 
tween a = 0.0° and 2.05° in which the flow is at- 
tached and both surfaces have positive pressure co- 
efficients. (See fig. 44.) The lift curve of these thick 
wings was found to be quite different from that of the 
flat wings. In particular, the linearity of the lift curve 
for all wings was found to be limited to a few degrees 
angle of attack. Moreover, the data of figure 41 show 
that the bounds established for the development of 
nonlinear lift for thin flat wings are not applicable 
to thick wings because significant amounts of non- 
linear lift are evident well within the flat-wing linear 
lift region. An examination of the surface pressure 
coefficient data indicated that the increased lift ef- 
fect may be attributed to the wing upper surface. 
At zero lift the thick-wing upper surface has large 
positive pressures, and as we have already seen from 
the flat-wing data, there is a minimum upper sur- 
face pressure which may be attained at a given Mach 
number. The differences between these two pressure 
levels define the maximum attainable upper surface 
wing loading. If we compare this value to that for 
a zero-thick or flat wing, at zero lift the flat wing 
has a pressure coefficient value of approximately 0.0 
compared with a value of 0.05 to 0.10 for the thick 
wing. (See fig. 44.) This difference would reduce the 
allowable change in upper surface pressure coefficient 
for the flat wing and, thus, the resultant wing lifting 
potential. These effects are quantified in figures 42 


14 



through 44 for changes in thickness, Mach number, 
and angle of attack. 

The lower surface pressures show an equal effect 
of both the wing leading-edge angle and wing thick- 
ness. Mach number effects are shown in figure 43 for 
the 30° biconvex wing geometry at an angle of attack 
of 8° . The data show trends similar to those observed 
for the flat delta wings, a reduced upper surface load- 
ing with increased Mach number, and a general in- 
sensitivity to Mach number on the lower surface. A 
comparison of the thick- wing lower surface data with 
the flat-wing data (fig. 29) shows a reversed variation 
with increasing Mach number. The flat-wing data in- 
dicate a slight increase in lower surface pressure with 
increasing Mach number whereas the thick-wing data 
show a decrease in pressure. This pressure reduction 
is actually misleading, because the pressure data at 
an angle of attack of 0° also show a pressure reduc- 
tion with increased Mach number, which occurs at a 
rate equal to that observed at angle of attack; this 
results in an invariant lower surface loading. 

Presented in figure 44 is the effect of angle of 
attack on the spanwise pressures for the 30° biconvex 
geometry at M = 1.30. The data show trends 
similar to those for the flat wings; however, the levels 
are quite different especially at an angle of attack 
of 0°. As mentioned previously, the positive upper 
surface pressure coefficient at a = 0° for the thick 
wings allows for a greater increment in upper surface 
loading with increased angle of attack. In particular, 
the data show an average pressure coefficient of 0.075 
at an angle of attack of 0°, which is approximately 
9 percent of the maximum flat-wing upper surface 
lift force, with total vacuum pressure coefficient on 
the wing upper surface being assumed. 

The distribution of loading between the upper 
and lower surfaces of the 30° biconvex wing at a = 8° 
and 16° for Mach numbers of 1.30, 2.00, and 2.80 is 
shown in figure 45. The data show a gradual shifting 
in lift from the upper surface to the lower surface 
with increasing Mach number and angle of attack. If 
the thick-wing results of figure 45 are compared with 
the flat- wing data of figures 29 through 31, we see a 
delay of the dominance of the lower surface on the 
total lift with increasing Mach number. 

A summary of the individual normal-force char- 
acteristics of the upper and lower surfaces for the 
three thick delta wings is presented in figures 46, 
47, and 48. Figures 46 and 48 were formulated in 
the same manner as those for the flat wings. Thick- 
wing upper surface normal-force data are presented 
in figure 46 along with the flat-wing curves from 
figure 34 for Mach numbers 1.30, 2.00, and 2.80. 
The thick-wing data collapse into three groups of 


constant Mach number in which the data agree 
reasonably well over the range of the correlation 
parameter &n(3 cot A. At low values of the correla- 
tion parameter, the thick-wing results fall below the 
flat-wing curves, and as the correlation parameter is 
increased the thick-wing data increase more rapidly 
and eventually cross over the flat-wing curves. The 
nonlinearity in the thick-wing data curves can be 
attributed to the large positive pressure coefficients 
which exist at a = 0°. Despite these differences be- 
tween the two data sets, some of the general trends 
of the thick wings agree reasonably well with the flat- 
wing results. 

The effect of the pressure at a = 0° is high- 
lighted in figure 47, which presents the upper sur- 
face minimum pressure coefficient at M = 1.30, 2.00, 
and 2.80. The data show a significant delay in the 
occurrence of a negative pressure coefficient due to 
increased thickness and leading-edge slope. This ef- 
fect is most evident for the data at M = 1.30 where 
there is a 0.25 variation in the minimum pressure co- 
efficient between the three wings at an angle of attack 
of 12°. 

The effect of thickness on the lower surface 
normal-force coefficient compares well to the flat- 
wing results of figure 36 (fig. 48); this indicates that 
lower surface normal-force characteristics are inde- 
pendent of wing thickness and Mach number. 

Camber effects. In designing a wing for efficient 
supersonic flight, the geometric parameter which is 
usually optimized after the wing planform selection 
is the camber and twist distribution. The purpose 
of this section is to evaluate the following: the 

wing lift-curve slope sensitivity to camber, the effect 
of aspect ratio on nonlinearities in the flow, and 
the lift distribution between the upper and lower 
surfaces relative to that for a flat wing. These points 
were addressed by analyzing a parametric set of six 
conical, zero-thickness, cambered delta wings tested 
at M = 1.90. The geometries consisted of 15° 
streamwise leading-edge deflections of the outboard 
10 percent and 20 percent of the semispan for each 
of the three leading-edge sweep angles of 75°, 67.5°, 
and 58.25° (ref. 24). These wings were used to 
evaluate the effect of camber only, since no twist 
was applied to the spanwise sections. The data are 
presented in a similar fashion as those for thickness 
effects in figures 49 through 57. 

Unlike the thick-wing data which showed a uni- 
form progression in the location of the vortex with 
all parameters, the cambered wing data were very 
erratic (fig. 49). A review of the pressure data re- 
vealed that the lee-side flow characteristics were quite 
complex due to the sudden expansion about camber 
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hinge line. (See figs. 52, 53, and 54.) However, 
once flow separation occurs at the wing leading edge 
(a > 10°), the flow behaves similar to that for the 
thin flat wings. The data also show that increasing 
wing camber delays the formation of a wing leading- 
edge vortex. These results are analogous to the effect 
of increasing the leading-edge sweep angle for thick 
wings. However despite these flow complexities the 
agreement between the flat-wing and cambered-wing 
lift-curve slopes was quite surprising (fig. 50). The 
cambered-wing data show a variation in lift-curve 
slope similar to that observed for thick wings with an 
increased lift-curve slope for /? cot A greater than 0.5. 
This increased lifting efficiency at the higher values 
of /? cot A may be attributed to the increased load- 
ing potential of the wing leading edge. The cambered 
wing at an angle of attack of 0° would have a large 
negative load on the deflected portion of the wing 
leading edge. The magnitude of this load at a = 0° 
would be equivalent to the additional lift increment 
available to the cambered wing compared with the 
flat wing. This effect is supported by the compar- 
ison of flat and cambered nonlinear lift increments 
which show the existence of nonlinear lift for all val- 
ues of the correlation parameter (fig. 51). To further 
evaluate the local aerodynamics of cambered delta 
wings, detailed upper and lower surface pressure co- 
efficient (at x/c r ss 0.66) distributions highlighting 
the individual effects of camber, leading-edge sweep, 
and angle of attack are presented in figures 52, 53, 
and 54, respectively. 

The effect of camber on the spanwise surface 
pressure distributions is presented in figure 52 for 
the 75° swept wing at an angle of attack of 8°. 
A significant variation in the leading-edge upper 
surface pressures is evident due to increasing camber; 
however, the loadings are equivalent. On the lower 
surface the spanwise pressure distributions show less 
sensitivity to camber. 

Presented in figure 53 is the effect of leading-edge 
sweep on the spanwise pressure distributions at an 
angle of attack of 8°. The data show that the major 
influence of leading-edge sweep is an increase in the 
lower surface loading with decreasing sweep; this also 
correlates well with the flat-wing data of figure 32(b). 
There is also a significant influence on the position 
of the vortex. 

The variation in the spanwise surface pressure 
distributions with angle of attack for cambered wings 
is presented in figure 54. The data clearly show the 
loading at the wing leading edge for an angle of attack 
of 0°, which was alluded to earlier. The integrated 
lift increment between the upper and lower surface 
spanwise pressure distributions at a = 0° provides 
an additional increment in lift to the wing similar to 


that observed for the thick wings. Also evident are 
the large hinge-line-induced expansion pressures on 
the wing upper surface (a = 8°) and the oscillatory 
nature of the lower surface pressure due to the irreg- 
ular geometry. As with the thick wings, an integra- 
tion of the pressures to obtain the upper and lower 
surface wing loadings resulted in values comparable 
with those for the flat wings (figs. 55 and 57). 

The upper surface normal-force data for all cam- 
bered wings are presented in figure 55. The data are 
shown to collect near or below the flat-wing data at 
M = 1.90 from figure 34 (solid line) and show little 
effect of leading-edge sweep or camber. Presented 
in figure 56 are the upper surface minimum pressure 
coefficients for all cambered wings. The data fall into 
a broad band which lies below and parallel to the flat- 
wing data. This large scatter of the CJ} min data can 
be attributed to the interaction between the hinge- 
line expansion and the leading-edge separation and 
produces both a wide range and a dramatic change 
in the lee-side flow with changes in camber, wing 
sweep, and angle of attack. The effect of camber 
on the lower surface normal-force coefficient is pre- 
sented in figure 57. The lower surface characteristics 
are shown to separate according to the wing leading- 
edge sweep as indicated by the flat-wing data (solid 
line), and the data show little effect of camber. A 
comparison of the two data sets seems to indicate a 
reduction in the lower surface normal force due to 
camber; however, the cambered delta wings of refer- 
ence 24 were constructed by a deflection of the out- 
board segment of the reference flat wing. This would 
result in leading-edge sweep angles for the cambered 
wings which are greater than those listed. 

Comments on lifting characteristics. The aero- 
dynamic characteristics of delta wings at lifting con- 
ditions have been evaluated for the effects of wing 
leading-edge sweep, leading-edge bluntness, and wing 
thickness and camber and then summarized in the 
form of graphs which may be used to represent the 
aerodynamic characteristics of delta wings. Empiri- 
cal correlation curves derived from experimental data 
have been developed for the lift-curve slopes, nonlin- 
ear lift effects, maximum lift, longitudinal stability, 
and distribution of lift between the upper and lower 
surfaces of a wing. However, the impact of airfoil 
thickness, maximum thickness position, leading-edge 
bluntness, wing leading-edge sweep, and lift coeffi- 
cient on the drag and the drag-due-to-lift character- 
istics is shown theoretically. 

Real-Flow Wing Design 

Wing design studies at supersonic speeds have 
typically been prefaced with the terms “linear” or 
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“nonlinear” depending upon the theoretical method 
which was employed within the investigation. In ad- 
dition, the previous nonlinear studies have employed 
traditional linear-theory rules in planform and airfoil 
selection. Previous supersonic wing design studies 
have also typically employed airfoils with thicknesses 
between 3 and 6 percent and the maximum thick- 
ness located at 40 percent to 50 percent chord. The 
aerodynamic reasoning for selecting or limiting the 
selection of wing geometries to these geometries has 
never been documented. In fact, linear- theory aero- 
dynamic graphs which are typically used in prelimi- 
nary design indicate that significant drag reductions 
may be attained with alternate airfoil shapes. (See 
fig. 6.) The intent of this section of the paper is to 
employ the empirically derived graphs of the previ- 
ous sections in a systematic fashion to select wing 
geometric characteristics conducive to high levels of 
aerodynamic performance. 

The zero-lift, low-lift, and high-lift aerodynamic 
characteristics previously presented for delta wings 
show a large influence due to changes in wing ge- 
ometry and flow conditions. At zero lift, theoreti- 
cal analysis with an inviscid nonlinear method shows 
significant variations from the standard linear-theory 
curves due to changes in airfoil shape, airfoil thick- 
ness, leading-edge bluntness, and wing leading-edge 
sweep. A review of the zero-lift-curve drag analysis 
shows that for a cruise or low-lift-dominated design 
the diamond airfoil could provide 20-percent to 50- 
percent reductions in wave drag compared with the 
NACA modified four-digit-series airfoil. However, if 
the wing is to operate over a range of Mach number 
and lift coefficient the problem becomes more com- 
plex and a blunt airfoil should be employed to im- 
prove the drag-due-to-lift characteristics. The zero- 
lift wave-drag analysis for the NACA modified four- 
digit-series airfoil (fig. 13) shows that for values of 
/? cot A between 0.6 and 0.8 locating the airfoil max- 
imum thickness at the 20-percent chord provides low 
drag, and the analysis presented in figure 26 shows 
an improvement in the drag-due-to-lift characteris- 
tics for the same geometry. The design of a wing 
within this range of / 3 cot A would provide the op- 
portunity to achieve good aerodynamic performance. 

A review of the lift-curve-slope data of figure 20 
and the computed drag-due-to-lift results of figure 27 
supports the range of (3 cot A from 0.6 to 0.8 as a fea- 
sible design region. The data of figure 20 show that 
at a value of /? cot A of 0.6 the data for the sharp 
leading-edge wings level off and at a value of /? cot A 
of 0.8 the same occurs for the blunt leading-edge 
data. For flat wings, the lift-curve slope is inversely 
proportional to the drag-due-to-lift characteristics; 
therefore, a leveling off of the lift efficiency should 


indicate the point of optimum performance. How- 
ever, the data of figure 27 show that the drag-due-to- 
lift characteristics are much more complex. Whereas 
the data of figure 20 showed no distinguishable effects 
of aspect ratio, the drag data show large variations 
in performance due to both wing sweep and lift co- 
efficient. The data show that the value of /? cot A 
at which the drag-due-to-lift characteristics reach a 
minimum varies between 0.4 and 0.8 depending upon 
lift coefficient and aspect ratio. A review of the zero- 
lift wave-drag data (fig. 13) shows that a value of 
( 3 cot A of 0.4 corresponds to the bucket of the curve 
for m = 0.2; as a result, the design space is ex- 
panded to values of /? cot A between 0.4 and 0.8 for 
completeness. 

The typical application of thickness to uncam- 
bered delta wing results in a geometry which is con- 
ical about the wing tip. This classical application 
of thickness to a swept wing is no doubt less than 
optimum for supersonic flight. Experimental data 
(ref. 38) and theoretical analysis (ref. 2) show that 
the flow over a swept wing at supersonic speeds tends 
to be conical about the wing apex and not coni- 
cal about the wing tip as observed in subsonic flow 
for unswept wings. The conical nature of the flow 
field over the wing produces favorable and unfavor- 
able pressure fields on the wing surface. For a wing 
at lift, the flow over the wing upper surface would 
be characterized by an expansion over the leading 
edge which is followed by a recompression to a more 
positive pressure as the flow moves inboard. The 
location of the recompression has been observed to 
lie along a ray emanating from the wing apex, in- 
dependent of the wing geometry. If the upper sur- 
face is divided into four quadrants, defined by the 
intersection of the airfoil maximum thickness line 
and the recompression line, two favorable and two 
unfavorable performance regions may be identified. 
The two unfavorable regions, which contribute to the 
drag, are the inboard forward portion and the out- 
board aft portion of the wing. The inboard forward 
portion of the wing experiences a recompression of 
the flow prior to the airfoil maximum thickness line; 
this results in more positive pressures acting on a 
forward-facing wing surface. On the other hand, the 
outboard aft portion of the wing is characterized by 
a rearward-sloping surface which combines with the 
high negative pressure coefficients to produce high 
drag levels. The other two quadrants of the wing up- 
per surface would have pressure fields which combine 
favorably with the local surface geometry to produce 
drag reductions. These observations suggest that im- 
proved supersonic performance could be achieved if 
the wing designer configures the wing geometry to 
take advantage of the natural conical flow structure 
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against it. Future supersonic wing designs should 
break from tradition and begin to explore new airfoil 
shapes and distributions of airfoils, such as increas- 
ing leading-edge blunt ness, airfoil maximum thick- 
ness, and airfoil maximum thickness position with 
increasing spanwise position to take full advantage 
of the conical nature of the flow. 

To define a high-lift (0.2 < Cl < 0.4) wing de- 
sign philosophy, an understanding of the flow lim- 
itations which exist is required. At these high-lift 
conditions, the flow will undoubtedly separate at the 
leading edge for an uncambered subsonic leading- 
edge wing independent of airfoil shape. The lee- 
side flow condition which would exist for delta wings 
at angle of attack is shown in figure 17. The data 
of figure 17 clearly show that within the feasible 
range of values of / 3 cot A identified for low lift 
(0.4 < (3 cot A < 0.8), Mjv < 1.00, the lee-side flow 
is characterized by a leading-edge separation. By 
limiting the high-lift design to separated flows and 
by imposing the minimum supersonic Mach number 
identified in figure 28 and an arbitrary selection of 
A = 75° as the maximum sweep angle, the M N and 
aw space can be reduced significantly. (See fig. 58.) 
The primary premise on which the high-lift wing de- 
sign concept is to be based is the assumption that 
the lift force on the wing upper surface must be equal 
to or greater than the lift force on the lower surface. 
The supporting argument for these criteria is that the 
negative pressure coefficients acting on the upper sur- 
face at the leading edge of the wing are the primary 
mechanism for creating aerodynamic thrust; thus, if 
the upper surface of the wing does not produce at 
least 50 percent of the lift, then significant drag re- 
ductions may not be realized. In order to satisfy 
these conditions, delta wing upper surface and lower 
surface normal-force coefficient design graphs have 
been developed from the data of figures 34 and 36 
(figs. 59 and 60). The upper surface normal-force 
coefficient characteristics are strongly dependent on 
Mach number; for example, an increase in Mach 
number from 1.50 to 2.00 reduces the upper surface 
normal-force coefficient by 50 percent. In a similar 
manner the lower surface normal-force coefficient is 
a strong function of leading-edge sweep. Additional 
data show that both the upper surface and lower sur- 
face normal-force coefficients of delta wings are not 
significantly influenced by thickness or camber; thus, 
the application of these graphs to the development of 
a high-lift wing design space for delta wing geome- 
tries should be acceptable. 

By imposing a design Cjy value and requiring that 
the design Cjy can be distributed between the upper 
and lower wing surfaces such that the upper surface 
normal-force coefficient always equals or exceeds the 


lower surface normal-force coefficient (Cjy > C^ N ), 
the data presented in figures 59 and 60 can be 
used in an iterative sense to define a feasible high- 
lift separated-flow design space. To assist in the 
iteration process, plots of Afjy and a^y against a are 
presented in figure 61. Presented in figure 62 is a 
matrix of points, identified by their associated Mach 
number and leading-edge sweep values, that satisfy 
the distributive normal-force coefficient relationship 
(Cjy > Cfl) for a design condition of Cjy « 0.40. 
Each of these points in the matrix was determined 
by iterating through figures 59 and 60 until the 
distributive requirement was satisfied. A sample 
design iteration for C ^ = 0.4 and A = 75° is initiated 
by determining the o/y at which = 0.2 and 
then through trial and error the maximum Mach 
number is determined at which Cjy = 0.2. Once 
aw an d M are known, the graphs in figure 61 can 
be used to compute A4)y for plotting in figure 62. 
The matrix of feasible solutions presented in figure 62 
can be thought of in terms of lines of constant Mach 
number (solid lines) in which both a and A vary and 
lines of constant leading-edge sweep (dashed lines) 
in which only M varies. The boundaries of a given 
design space are defined on the right by the maximum 
wing leading-edge sweep under consideration, on the 
bottom by the minimum Mach number, and on the 
left by the minimum allowable wing leading-edge 
sweep under consideration. The upper boundary is 
defined by iterating through the design process to 
determine the maximum Mach number for which a 
given wing geometry (leading-edge sweep) satisfies 
the distributive requirement. Shown in figure 62 are 
the maximum Mach numbers for several leading-edge 
sweep conditions. 

Presented in figure 63 is the design space for de- 
sign values of Cj y of 0.40 and 0.20. The data indicate 
that a reduction in the design value of C jy expands 
the range of feasible solutions and shifts it to a lower 
value of Gfjv and Mj y . A decrease in design C \ y allows 
for the extension to higher Mach numbers. The ap- 
parent reduction in the size of the design space with 
reduced design Cj y is strictly graphical in nature due 
to the nonlinearly decreasing relationship between a 
and ajsf. The design-space concept for moderate- 
to-high-lift conditions (0.2 < Ci < 0.4) was evalu- 
ated with data from previous supersonic wing designs 
(ref. 39). The study showed that designs conducted 
for conditions which lie within the appropriate design 
space had lower drag due to lift compared with those 
outside the design space. 

To further evaluate the design-space concept and 
to relate the high-lift requirements to the low-lift 
requirements, lines of constant values of (3 cot A 
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of 0.4, 0.6, and 0.8 are presented in each plot. The 
value of for each point in figure 63 corresponds to 
the value of a jy at which C^ N equals 50 percent of the 
total design C jy value for each wing sweep. The value 
of Mjv was then determined by iterating through fig- 
ure 61(a) for the particular value of /? cot A. The 
graphs show that decreasing both f3 cot A and the 
design C n expands the design matrix. For the de- 
sign (7/v of 0.4, decreasing (3 cot A from 0.8 to 0.6 
and 0.4 increases the maximum wing leading-edge 
sweep for designs from 55° to 65° and to a value 
greater than 75°, respectively. Similarly decreas- 
ing design Cjy to 0.2 increases the maximum wing 
leading-edge sweep from 55° to 60° for a design at 
a value of f3 cot A of 0.8. These results have been 
summarized in figure 64 in which the limiting sweep 
angle for a delta wing has been plotted against the 
leading-edge sweep parameter j3 cot A for design val- 
ues of Cjsj of 0.2 and 0.4. This graph clearly shows 
the dominance that the (3 cot A condition has on the 
ability to produce efficient lifting conditions required 
for high-lift wing design. 

If the high-lift design trade data of figure 64 are 
combined with the zero-lift data of figure 13 and the 
low-lift data of figures 20 and 27, it may be con- 
cluded that a value of /3 cot A of approximately 0.6 
composed of A = 65° and M — 1.63 would pro- 
vide the optimum performance for lift coefficients 
between 0.0 and 0.4. For a conventional design, 
the wing would be configured with a moderately 
thick airfoil ( t/c « 0.04) with maximum thickness 
located at 20 percent of the chord to minimize the 
zero-lift drag at (3 cot A = 0.6 and to maximize the 
low-lift drag-due-to-lift characteristics. (See fig. 26.) 
The 65° swept delta wing with aspect ratio of 1.86 
should provide a minimum zero-lift drag penalty and 
a significant improvement in drag due to lift com- 
pared with a more slender geometry. At M — 1.63, 
the effect of vacuum limit is minimal, providing a 
70-percent increase in upper surface lifting potential 
over the 75° design (/? cot A = 0.45) and an 8-percent 
decrease from the 55° design (/ 3 cot A = 0.8). If a 
simple variable camber system (simple flap), which 
is conical in nature, is added to the described de- 
sign, the aerodynamic performance across the range 
of lift coefficient could approach that of the practical 
goal established in reference 40. (See fig. 65.) Also 
presented in the figure are the flat-wing aerodynamic 
characteristics for delta wings with A = 1 and 2 from 
figure 27. The projected performance levels proba- 
bly represent the minimum and not the maximum 
allowable performance levels. If both the airfoil and 


the camber are integrated together to take full ad- 
vantage of the conical flow characteristics over the 
wing at zero-lift to high-lift conditions, then aerody- 
namic performance levels which exceed the practical 
goal limits might be realized. 

Concluding Remarks 

Through the empirical correlation of experimen- 
tal data and theoretical analysis, a set of graphs have 
been developed which quantify the inviscid aerody- 
namics of delta wings at supersonic speeds. 

The zero-lift wave-drag characteristics of delta 
wings with diamond, circular-arc, and NACA modi- 
fied four-digit-series airfoils were determined through 
the application of a nonlinear computational tech- 
nique. The nonlinear analysis varied substantially 
from the exact linear-theory predictions for all com- 
binations of geometry and flow parameters under 
study. The nonlinear analysis showed that for slender 
wings with highly subsonic leading edges the zero- 
lift wave-drag correlation relationship is maintained; 
however, as the wing geometry becomes nonslender, 
the flow about the wing becomes nonlinear and the 
relationships which define the zero-lift wave-drag cor- 
relation parameter are not maintained. 

The aerodynamic characteristics of delta wings at 
lifting conditions have been evaluated for the effects 
of wing leading-edge sweep, leading-edge bluntness, 
and wing thickness and camber and then summa- 
rized in the form of graphs which may be used to 
assess for the aerodynamics in the preliminary de- 
sign process. Empirical curves have been developed 
for the lift-curve slope, nonlinear lift effects, maxi- 
mum lift, longitudinal stability, and distribution of 
lift between the upper and lower surfaces of a wing. 
In addition, the impact of various airfoil parame- 
ters, wing leading-edge sweep, and lift coefficient on 
the drag-due-to-lift characteristics has been shown 
theoretically. 

The various graphs which detail the aerodynam- 
ics of delta wings at both zero-lift and lifting con- 
ditions were then employed to define a preliminary 
wing design approach in which both the low-lift and 
high-lift design criteria were combined to define a 
feasible design space. 


NASA Langley Research Center 
Hampton, Virginia 23665-5225 
December 11, 1987 
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Appendix 

Nonlinear Methodology 

The nonlinear computational technique (ref. 2) selected for analysis solves the nonconser- 
vative finite-difference analog of the full-potential equation in a spherical coordinate system. 
The method marches in the radial direction to obtain three-dimensional cross-flow solutions of 
the wing geometry. The ability to compute the flow over a wing with this method is restricted 
to conditions in which the flow in the marching direction remains supersonic; thus, the code 
would be expected to be more successful for higher free-stream Mach numbers. The solution 
technique is graphically depicted in figure Al. A series of three transformations are used in the 
code: the geometry in Cartesian coordinates is transformed to spherical coordinates followed 
by a stereographic projection (not shown) and a conformal mapping to a circle. A shearing 
transformation is then performed to the computational plane. Grid points are internally com- 
puted within the code and are positioned between the body surface (inner boundary) and the 
bow shock (outer boundary). The method has been shown to be ideally suited for the solution 
of three-dimensional wing flow fields (ref. 2). 

Solutions were obtained by fitting the bow shock and capturing all internal shocks. Due to 
the supersonic nature of the wing trailing edge, the wake was modeled as an extension of the 
wing surface. Sensitivity studies were conducted to determine marching step size and cross- 
plane grid resolution. These studies indicated a dependence of both parameters to free-stream 
flow conditions as well as geometry. Selection of the marching step size was based upon the 
analysis of a series of airfoil geometries at various flow conditions. A representative sampling 
of this analysis is presented in figure A2 in which the axially integrated drag values for three 
airfoil types are presented for step sizes of 0.2, 0.5, and 1.0, which correspond to 50, 20, and 10 
marching steps, respectively. Results of the analysis show a convergence of the solution with 
decreasing step size. Based upon these results and considerations of computational cost, a step 
size of 0.5 was selected for the study. Cross-plane mesh density was selected to be 30 x 30 based 
upon a similar parametric study. Presented in figure A3 are representative 30 x 30 cross-plane 
computational grids for the four basic airfoil types under investigation: diamond, circular arc, 
sharp NACA modified four-digit series, and blunt NACA modified four-digit series. The grid 
plots clearly show clustering of the grid points at the wing leading edge in both the radial and 
circumferential directions; this is critical to resolving the wing flow field. 

NCOREL A = 1 .0 delta wing 
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Table I. Experimental Data Set 



Planform 

Airfoil 

Flow conditions 

Data type 

Reference 

A 

Type 

f/c 

LE 

M 

a, deg 

NrJc 

Force 

Pressure 

Flow visualization 

3 

3 wings 
76° to 84.6° 

Modified 

biconvex 

4% at 
0.59c 

Sharp 

1.96 to 
3.30 

0 to 
47 

1.68 x 10® to 
4.48 x 10 6 

v/ 



4 

11 wings 
45° to 80° 

Diamond 

8% at 
0.18c 

Sharp and 
blunt 

1.62 to 
2.40 

—6 to 
6 

0.46 x 10 6 to 
1.39 x 10 6 

V 




4 wings 
65° to 75° 

Slab 


Sharp and 
blunt 

1.92 

—6 to 
6 

0.46 x 10 6 to 
1.39 x 10 6 

v/ 



5 

4 wings 
48° to 73° 

Diamond 

6% at 
0.20c to 
0.55c 

Sharp 

1.62 to 
2.41 

-2 to 
6 

2.36 x 10® to 
4.32 x 10 6 




6 

3 wings 
63.43° to 
82.87° 

Slab 


Sharp 

1.60 to 

4.60 

-5 to 
60 

1.41 x 10® to 
2.82 x 10 6 




7 

2 wings 
60° 

Diamond 

5% 

Sharp 

4.07 

-4 to 
50 

5.3 x 10® 

V 




71.6° 

Diamond 

8% 

Sharp 

4.07 

-4 to 
50 

6 x 10® 




8 

1 wing 
68.6° 

NACA-00 

4 to 
6.24% 

Blunt 

2.41 

0 to 
20 

1.04 x 10® to 
18.3 x 10® 

v/ 

v/ 


9 

3 wings 
45° and 64° 

Slab 


Sharp 

1.55 to 
2.32 

0to 

50 

0.18 x 10® to 
0.49 x 10® 

v/ 




Cn 

00 

0 

Diamond 

10% at 
0.50c 

Sharp 

2.40 




V 


10 

2 wings 
65° 

RAE 101 
Biconvex 

4% 

Blunt 

1.51 to 
2.0 

—6 to 
10 

3 x 10® 




11 

8 wings 
36.4° to 
70.6° 

Slab 


Blunt 

1.43 and 
1.71 

-5 to 
5 

0.24 x 10® to 
0.62 x 10® 

v/ 



12 

2 wings 
63.43° 

Diamond 

5% at 
0.20c and 
0.50c 

Sharp and 
blunt 

1.53 

-2 to 
16 

0.75 x 10® 

V 



13 

2 wings 
60° 

Diamond 

8% at 
0.18c and 
0.60c 

Sharp 

4.04 

-15 to 
15 

6 x 10® 1 

V 



14 

2 wings 
68.4° 

NACA-00 

4 to 

6.24% at ! 
0.30c 

Blunt 

1.6 and 
1.9 

-4 to 
16 

2.2 x 10® to 
18.4 x 10® 

V 

V 


15 

1 wing 
63.26° 

NACA-0003 

6% 

Blunt 

1.30 to 
2.22 

-6 to 
18 

3.68 x 10® 

s/ 



16 

2 wings 
45° and 
63.4° 

Biconvex 

5% at 
0.50c 

Sharp 

1.45 and 
1.97 

0 to 
50 

0.68 x 10® to 
3.27 x 10® 




17 

3 wings 
45° to 
63.26° 

NACA-0003 
and Biconvex 

3 to 8% 

Sharp and 
blunt 

1.2 to 
1.7 

-4 to 
20 

1.4 x 10® to 
3.0 x 10® 

V 
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Table I. Concluded 


Reference 

Planform 

Airfoil 

Flow conditions 

Data type 

A 

Type 

t/c 

LE 

M 

a, deg 

NRe/c 

Force 

Pressure 

Flow visualization 

18 

3 wings 
53° to 
66.6° 

53° 

65A003 

Flat 

3% 


1.62 to 
2.41 

1.94 and 
2.41 


2.3 x 10 e to 
4.65 x 10 6 

7 

v/ 


19 

3 wings 
63.43° to 
82.87° 

Slab 


Sharp 


0 to 
55 


■ 




6 wings 
58.25° to 
85° 

Flat 


Sharp 

1.90 


S3QnQQfijQ|| 


V 



4 wings 
52.5° to 
75° 

Flat 


Sharp 



0.86 x 10 e to 
2.46 x 10 6 


V 


22 

2 wings 
45° 

45° 


6% 

Blunt 

Sharp 



1.8 x 10 6 


V 


23 

1 

Cross sections: 
30° Diamond 
30° Biconvex 
60° Biconvex 

Conical 

Sharp 

1.30 to 
2.80 

0 to 
16 

3.63 x 10 B to 
2.26 x 10 6 


V 


24 

3 wings 
58.25° to 
75° 

Flat 

Cambered 

Sharp 

1.90 


H| 


sj 


25 

8 wings 
71.56° 

Slab 

Cambered 

Sharp 

1.40 and 
1.80 

~8 to 
12 

0.66 x 10° 

V 
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Figure 6. Linear- and nonlinear-theory-predicted zero-lift wave drag for delta wings with diamond airfoils. 




Figure 7. Spanwise pressure distributions predicted by nonlinear theory for delta wing with A = 1.0 and 
diamond airfoil with r = 0.04 at M = 2.60. 
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Figure 10. Linear- and nonlinear-theory-predicted zero-lift wave drag for delta wings. 
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Figure 13. Nonlinear-theory-predicted zero-lift wave drag of delta wing with NACA modified four-digit-series 
airfoils. 
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Figure 20. Nonlinear-theory-predicted lift characteristics of thin delta wings. 


- 20 Ci 

20° a o° 

Airfoil 


Blunt 

Sharp 

Reference 


□ 

6 


> 

32 


D 

14 

♦ 


15 


O 

16 

W 


17 


0 

18 


o 

19 





46 



m | a = 20 



P cot A 

Figure 24. Nonlinear pitching-moment characteristics of thin delta wings 
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Figure 25. Wing panel center-of-pressure location for thin delta wings (ref. 6). 
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Figure 26. Nonlinear- and linear-theory-predicted drag-due-to-lift characteristics of delta wings with A = 1.0 
and modified four-digit-series airfoils at M = 1.41. 
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(c) M = 1.20 (ref. 19). 

Figure 28. Continued. 
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(e) M = 2.16 (ref. 6). 
Figure 28. Continued. 
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Figure 30. Effect of angle of attack on spanwise surface pressure distributions of delta wing with A 
for z/c r = 0.625. 
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Figure 33. Effect of M and A on the distribution of lift for delta wings. 
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(a) Reference 6. 

34. Wing upper surface normal-force coefficient for thin delta wings. 
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(c) Reference 20. 
Figure 34. Concluded. 
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(b) Reference 21. 
Figure 35. Continued. 


M = 1.90 



0 10 20 30 40 50 


a |sj P cot A 

(c) Reference 20. 
Figure 35. Concluded. 
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Figure 38. Effect of iV/j e on C“ min for delta wings with A = 68.6°. 
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Figure 39. Location of vortex action line for thick delta wings (ref. 23). 
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Figure 40, Lifting characteristics of thick and thin delta wings (ref. 23). 







Figure 42. Effect of thickness on spanwise pressure distributions of thick delta wings at M 
and ol = 8.0° (ref. 23). 
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Figure 44. Effect of angle of attack on spanwise surface pressure distributions of thick delta wings 
at M = 1.30 (ref. 23). 
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Figure 46. Comparison of upper surface normal-force characteristics on thin and thick delta wings (ref. 23) 
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Figure 47. Comparison of upper surface minimum pressure coefficient on thin and thick delta wings (ref. 23). 
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Figure 49. Location of vortex action line for cambered delta wings (ref. 24). 
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Figure 51. Nonlinear lift effects of cambered delta wings at M — 1.90 (ref. 24). 








Figure 53. Effect of A on spanwise surface pressure distributions for cambered delta wings at M 
and a = 8° (ref. 24). 
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Figure 55. Upper surface normal-force coefficient for cambered delta wings (ref. 24) 
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Figure 59. Delta wing upper surface normal-force coefficient design curves. 


96 








2 2 1/2 
= M cos A (1 + sin a + tan A) 


M 

10.00 5.00 4.00 3.00 2.00 1.50 1.20 




« N 


Figure 61. Variation of Mjy and with a. 
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Figure 62. Wing design space concept. 



Figure 64. Identification of limiting leading-edge sweep angle. 
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Figure 65. Identification of practical aerodynamic performance goals. 
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